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Abstract 

We derive asymptotics for the quenched probability that a critical branching Brownian motion 
killed at a small rate e in Poissonian obstacles exits a large domain. Results are formulated in terms 
of the solution to a semilinear partial differential equation with singular boundary conditions. The 
proofs depend on a quenched homogenization theorem for branching Brownian motion among mild 
obstacles. 
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1 Introduction 



In the present work, we are interested in the long-term behaviour of branching Brownian motion killed 
in Poissonian obstacles. Let us start by describing a simple special case of our results. We consider 
a critical branching Brownian motion in M, d (d > 1), where all initial particles start from the origin. 
£T) ■ We assume that particles are killed at a (small) rate e > within random balls of fixed radius, whose 
, centers are distributed according to a homogeneous Poisson point process on M. d . Then, how many 
initial particles do we need so that, with high probability, one of their descendants reaches distance R 
from the origin ? Let p £ (R) be the (quenched) probability for our randomly killed branching Brownian 
motion starting with a single particle at to visit the complement of a large ball of radius R centered 
at the origin. The preceding question is equivalent to determining the limiting behaviour of p £ (R) 
^ ■ when e tends to and simultaneously R tends to infinity. 

The answer involves several regimes depending on the respective values of e and R. If e is small 
in comparison with 1/R 2 , the killing phenomenon does not matter and the result is the same as if 
there were no killing: p e (R) behaves like a constant times 1/R 2 (informally, the branching process 
must survive up to a time of order R 2 so that at least one of the particles travels a distance R, and 
well-known estimates for critical branching processes then lead to the correct asymptotics). On the 
other hand, if e is large in comparison with 1/R 2 , then the probability p £ {R) decreases exponentially 
fast as a function of Ryfe: See Proposition [1] below. 

Our main results focus on the critical regime where eR 2 converges to a constant a > 0. We show that 
the probability p £ {R) behaves like R~ 2 , as in the case without killing, but with a multiplicative constant 
which depends on a and can be identified as the value at the origin of the solution of a semilinear partial 
differential equation with singular boundary conditions. A key tool to derive these asymptotics is a 
quenched homogenization theorem which shows that our branching Brownian motions among obstacles, 
suitably rescaled, are close to super-Brownian motion killed at a certain rate depending on a. 

Let us formulate our assumptions more precisely in order to state our results. First, let us define 
the collection of obstacles. We denote the set of all compact subsets of M. d by K. This set is equipped 
with the usual Hausdorff metric djj- Recall that (/C, dn) is a Polish space. For every r > 0, fC r denotes 
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the subset of K, which consists of all compact sets that are contained in the closed ball of radius r 
centered at the origin. Let O be a finite measure on /C, and assume that is supported on /C ro for 
some r > 0. Let 

be a Poisson point measure on R d x K, with intensity A^ ® 0, where A^ stands for Lebesgue measure 
on M. d . We assume that this point measure is defined on a probability space (fi, P) and we denote the 
generic element of fi by w. Our set of obstacles is then defined by 



r w = \J{x i + K i ), (i) 



iel 

where obviously Xi + Ki = {z = x% + y : y G K{\. Note that we use the notation T ro to emphasize 
that the set of obstacles depends on the variable w representing the environment. Let us also define a 
constant k by 

k = P(0 G = 1 - exp ( - J Q(dK) AdOfif)) • 

To avoid trivial cases, we assume that k > 0, or equivalently 0(Ad(iT) > 0) > 0. By translation 
invariance, we also have P(x G r ro ) = k for every a; G 

Let us now introduce the sequence of branching Brownian motions of interest. Given w G and 
a parameter e > 0, we consider a branching Brownian motion on M. d such that 

• each particle moves around in M. d according to the law of Brownian motion killed at rate s within 

r • 

- 1 - VJ 1 

• each particle branches at rate 1. During a branching event, the particle generates a random 
number of offspring, according to an offspring distribution v which has mean one and finite 
variance a 2 G (0, oo). 

This branching Brownian motion is denoted by Z m ' £ = (Z™' £ )t>o, where Z^ 7 ' 6 stands for the sum of 
the Dirac point masses at the particles alive at time t. The processes Z™' £ are defined on a probability 
space Q. For every finite point measure \x on M. d , we use the notation for the probability measure 
on under which each of the processes Z™ ,£ starts from /i. 

Let A be a bounded domain of class C 2 in M. d containing 0. We say that the branching Brownian 
motion Z m ' e hits A c if there exists t > such that Z™' e (A c ) > 0. We are interested in asymptotics for 
the quantity 

Fs (Z™' £ hits {RA) C ) 

when R — > oo and e — > 0. Here we use the obvious notation RA = {z = Ry : y G A}. 

Theorem 1. For every a > 0, let ur a \ = (uu\(x),x G A) be the unique nonnegative solution of the 
singular boundary value problem 

1 a 2 

-Au=— u +au m A, ^ 



L \dA 



+oo . 



Then, 

lim ( sup i? 2 P 5o (Z ro ' £ hits (RA) C ) - u (k£R2) (0) )=0, P(dw) a.s. 

R-^oo \ e > / 

Let us state a corollary of the theorem, which is motivated by the simple question we asked at 
the beginning of this introduction. For every a > 0, we denote by u° a ^ the solution of the boundary 

problem (|2|) when A is the open unit ball of M. d . 
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Corollary 1. For every e G (0,1), let n £ be a positive integer. Assume that en £ — > b as e — > 0, 
for some b > 0. Denote by i? ro ' e the maximal distance from the origin attained by a particle of the 
branching Brownian motion Z ro,£ . Then, P(dro) a.s., the law of ^R m,£ under F n ,s converges as 
e — > towards the probability measure tv^, on R + defined by 

7Tb([0,r]) = exp(-^-u° Kr) (0)) 

for every r > 0. 

In the setting of Theorem [TJ it is not hard to see that it( a ) (0) tends to as a — > oo (see Lemma 
[2] below) and thus this theorem does not give much information when R — > oo and e — > in such a 
way that eR 2 tends to oo. In that case, the next proposition provides an exponential decay, which 
contrasts with the preceding theorem. Since our bounds are clearly not optimal, we consider only the 
case when A is a ball. We denote the open ball of radius r centered at the origin by 5(0, r). In the 
general case we may apply the bounds (i) and (ii) of the proposition after replacing A by a ball B(0, r) 
such that B(0,r) D A or B(0,r) C A respectively. 

Proposition 1. (i) There exists a positive constant Cq = Cq(v) such that, for every w G fi, R > 1 
and e G [1/R 2 ,1], _ 

F So (Z m < e hits B(0,R) c ) > C s exp(-i?v / 2e). 

(ii) There exists a positive constant C\ = Ci(0) such that for every R > 1 and e G [1/R 2 , 1], 

P®P 5o (Z ro ' e hits B(0,R) c ) < exp(-Cii?\/e)- 

Consequently, we can find two positive constants C2 = C2(0) and C3 = C^(Q) such that P(dw) a.s., 
for every sufficiently large R and every e G [C^log log R) 2 /R 2 , 1], 

P 5o (Z ro ' £ hits B(0,R) c ) < exp(-C 3 RVe). 

Part (i) of the proposition is derived from an estimate about branching Brownian motion killed 
homogeneously at rate e, which explains why this bound holds for every w G fi and does not depend 
on the measure 0. The bounds in (ii) follow from an estimate for Brownian motion killed in mild 
obstacles and therefore do not depend on the offspring distribution v. The first assertion in (ii) may 
be compared to Proposition 5.2.8 in |Szn98] , 

Proposition [1] only gives rather crude estimates, and it would be of interest to obtain more precise 
information on the decay of the quenched probabilities P^Z™' 6 hits (RA) C ) in the case when eR 2 
tends to 00. This leads to large deviations problems in the spirit of the work of Sznitman [Szn98j, 
which we do not address here. 

A major ingredient of the proof of Theorem [T]is the following quenched homogenization result. We 
need to introduce a rescaled version of the process Z ro,£ . For every e > and every t > 0, let us define 
a random measure X™ ' £ on R d by setting, for any nonnegative measurable function ip on K , 

(XT' 5 , <p)=ej Zp? t (dx) v(e 1/2 x). (3) 

Here and later, the notation (/i, cp) stands for the integral of the function ip against the measure fx, 
whenever this integral makes sense. 

For every real x > 0, [x] denotes the integer part of x. 

Theorem 2. Except for a P '-negligible set of values of w , the law of (X^ u ' £ )t>o under ¥t £ -iig Q converges 
weakly as e — > 0, in the Skorokhod sense, to that of a super- Brownian motion with branching mechanism 
ip( K )(u) := \u 2 + ku started at 6q. 
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The definition of super-Brownian motion with branching mechanism ipr K \ is recalled in Section 2 
below. 

As a hint of why Theorem [2] should be true, notice that for a given realization of the obstacles, the 
probability that a single Brownian motion starting from and killed at rate e within T ro is still alive 
by time t > is given by 



E 



expj-EjT l rw (&)dfl}], (4) 

where £ denotes standard dimensional Brownian motion. Let us focus on the integral within the 
exponential in (J4j) . Averaging over the law of the obstacles and using Fubini's theorem, we obtain for 
each t > 



E 



E 



flr^Qds]] =ee\ f P[6 G T m ]ds 
Jo n l Jo 



£Kt. 



We can thus guess, and easily prove, that the rescaled Brownian motion (i/i £ e -i£, t > 0), which 
is killed at rate 1 within y/e converges to Brownian motion killed at homogeneous rate At as 
e — > 0. Theorem |2] shows that an analogous convergence indeed holds in our more general framework 
of branching Brownian motions, for any fixed w contained in a set of P-probability one. 

Let us briefly explain how Theorem [T] is derived from Theorem [2j Consider a sequence (e n ,R n ) 
such that R n — > oo and e n K\ converges to a positive constant a. By a simple scaling transforma- 
tion, the probability P So (Z m ' £n hits (R n A) c ) coincides with P S() (X m ' £n hits (b n A) c ), where b n = el/ 2 R n 
converges to \/a. We can then use Theorem [2] to investigate the asymptotic behaviour of the latter 
hitting probabilities. This limiting behaviour involves the corresponding hitting probabilities for super- 
Brownian motion, which are known to be related to solutions of semilinear partial differential equations 
from the work of Dynkin |Dy91| |Dy93| . One difficulty in implementing the preceding idea comes from 
the fact that the convergence in Theorem [2] is not strong enough to ensure that hitting probabilities for 
the processes {X™ ,£ )t>o converge to hitting probabilities for the limiting process. Much of the proof 
of Theorem [T] in Section 4 is devoted to a precise justification of this property (Lemma [7]). 

To complete this introduction, let us mention that branching Brownian motion and superprocesses 
among random obstacles have been studied recently in several papers, including Englander and den 
Hollander [EdH03] and Englander |Eng08]. These papers concentrate on the case of supercritical 
branching, in contrast with critical branching which is considered here. See also the survey |Eng07| , A 
homogenization theorem related to Theorem |2] has been proved in [Veb09j for super-Brownian motion 
among hard obstacles, in the case when the intensity of the obstacles grows to infinity but their 
diameters shrink to 0. There is a huge literature about Brownian motion and random walks among 
(hard or mild) obstacles, and the reader may look at the book of Sznitman |Szn98j for additional 
references. 

The rest of this paper is laid out as follows. In Section (2j we introduce the basic notation and 
objects, and state several results about hitting probabilities for spatial branching processes we shall 
need in the sequel. Theorem [2] and Proposition [1] are proved in Section [3] Theorem [1] and Corollary [1] 
are then derived in Section [4] 



2 Preliminaries 
2.1 Notation 

We denote the set of all finite measures on M. d by A4f(M. d ). This set is equipped with the weak topology. 
We write A4. p (M. d ) for the subset of A4f(M. d ) which consists of all finite point measures on W 1 . 

If E is a metric space, we denote the set of all bounded continuous functions on the space E by 
C(E) and we let ||/|| stand for the supremum norm of / G C(E). We write C 2 (M. d ) for the set of 
all twice continuously differentiable functions on and C 2 (M. d ) for that of all bounded functions in 
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C 2 (M. d ) whose first and second derivatives are also bounded. An index + added to this notation means 
that we require the functions to be nonnegative. We equip C 2 (M. d ) with the topology induced by the 
seminorms where for every R > 




If E is a Polish space, we let Dg[0, oo) be the set of all cadlag paths with values in E, equipped 
with the Skorokhod topology. 

If x G M. d and r > 0, B(x,r) denotes the open ball of radius r centered at x, and B(x,r) stands 
for the corresponding closed ball. More generally, the closure of a subset F of M. d is denoted by F. 
Lebesgue measure on M d is denoted by A^. 

Finally, the notation £ = (Ct)t>o will stand for a standard Brownian motion in M. d , which starts 
from x under the probability measure P x . It will also be convenient to use the notation £ ro ' £ for 
Brownian motion in M. d killed at rate e in the set As usual, the value of £ ro ' e after its killing time 
is a cemetery point A added to W d , and we agree that all functions vanish at A. 

2.2 Super- Brownian motion 

Let a > and set ip^(u) = ^-u 2 + au, for every u > (the offspring distribution u, and thus the 
parameter a > are fixed throughout this work). Super-Brownian motion with branching mechanism 
-0(a) is the continuous strong Markov process with values in 7^4/(1^), whose transition kernels (Qt)t>0 
are characterized as follows: For every g G C+(M. d ) and every n G A4j(M. d ), we have for every t > 

J Q t (fJ,,dfjf) exp(-(//,c?)) = exp(-(fj,,V t g}), (5) 

where the function Ut(x) = Vtg(x), t > 0, x G M d , is the unique nonnegative solution of the semilinear 
parabolic problem 

— = iAu-W)(«) in(0,oo)xM d , 
uo = g ■ 

Let Y = (Yf)t>o be a super-Brownian motion with branching mechanism V(a)' started at [i G 
.M/(M d ). Then, for every 5 G C\ (R d ) 

e -(Y u9 ) _ e -(Yo, 9 ) _ J (y s , A A g + e-^ds (6) 

is a martingale. It is well known that this martingale problem and the initial value fi characterize the 
law of Y. This is indeed an application of the classical "duality method" (see in particular Chapter 
4 in |EK 86j). The nonlinear semigroup g — > Vtg provides a deterministic dual to super-Brownian 
motion, and the duality argument then shows that if a measure-valued process started from /i satisfies 
the preceding martingale problem, the Laplace functional of its value at time t must be given by the 
right-hand side of (|5]). See Section 1.6 of [EthOOj for more details. 

2.3 Branching Brownian motion among random obstacles 

In view of our applications (and in particular because we want to refer to some results of [Ch91j). 
it will be convenient to give a more formal description of the branching Brownian motions that were 
already introduced in Section [T] above. Recall that our offspring distribution u is assumed to be critical 
and that var(i^) = a 2 G (0,oo). The probability generating function of v will be denoted by T. 
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Let T be a Galton- Watson tree with offspring distribution v (see e.g. [LG05J). As usual, we view 
T as a random finite subset of 

oo 

U :=U N n , 

n=0 

where N = {1, 2, 3, . . .}, and N° = {0}. If v = . . . , v n ) G £A{0}, tne parent of v is denoted by 
t) = (wi,... ,f n -i) and we also use the notation t> -< v' to mean that v' is a descendant of v distinct 
from v. Consider a collection (e v ,v G C/) of independent exponential random variables with parameter 
1, which is also independent of T ■ We define for every v G U its birth time a v and its death time (5 V 
recursively by setting a = and [3 = e , and for every v G U\{0}, 

cx v j3y , cv v -\- e v . 

Let us now construct the spatial motions. Fix a starting point x G R rf , and consider a collection 
(B v ,v G [/) of independent standard Brownian motions in M. d (started from 0), independent of T and 
of (e v ,v G U). For every v G ?7, define the historical path lv v = (<^\0 < t < f3 v ) associated with v 
in the following way. First ujf = x + Bf for < t < /3 . Then, if v G U\{0}, set = uif for a11 
< t < a v and 

^ = ui v +Bl av fora v <t<{3 v . 

A branching Brownian motion (without killing in obstacles) starting from 5 X is obtained by setting for 
every t > 0, 

Z t = Yl 6 "h 

where the notation v ~ t means that a v <t</3 v . 

In this formalism, it is now easy to introduce killing in obstacles. Consider yet another independent 
collection (j v ) V £i/ of independent exponential random variables with parameter 1, and define for every 
w G fi, e > 0, and for every v G U 

C £ ■= infuse [a v ,/3 v ) : £ drI T J^ r ) > e~ l lv ^, 

where inf = oo. By setting 

Z7' £ = ^2 Ht<C™' S and Cj' e =oo, for every v'Hv} ^ , 

we obtain a branching Brownian motion killed at rate e in the obstacle set starting from 5 X . 
An obvious extension of the preceding construction allows us to obtain branching Brownian motions 
starting from any point measure [i G Ai p (M. d ). 

We now recall a special case of the classical convergence of rescaled branching Brownian motions 
towards super-Brownian motion. For our applications, we state the case where particles are killed at a 
constant rate homogeneously over JR d (this case is obtained from the preceding construction of Z m,£ by 
replacing r ro by W 1 ). In the next two statements, for every e > 0, = (Z^)t>o denotes a branching 
Brownian motion with offspring distribution v, where particles are killed homogeneously over M. d at 
rate e. As previously Z^ starts from /i under the probability measure P^, for every /x G A4 p (B. d ). 

Proposition 2. Let a > 0. For every e > 0, define a measure-valued process (X\ ;)t>0 by setting 

{X l f\v) = e f Z^dx)^ 1 /^), 

for every tp G C(R d ). For every fixed i] > 0, the law of X^ under P^-ii^ converges as e — > 0, in the 
Skorokhod sense, towards the law of super-Brownian motion with branching mechanism ?/>( a ) starting 
from i]5q. 
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A proof of Proposition [2] can be found in Chapter 1 of [EthOOj in the case a = 0, and arguments 
are easily adapted to cover the general case. 

Finally, we shall use an estimate for the probability that a branching Brownian motion starting from 
5q exits a large ball centered at the origin. Similar estimates can be found in Sawyer and Fleischman 
|SF79| . but we provide a short proof for the sake of completeness. 

Lemma 1. Suppose that d = 1. There exist two positive constants Cq = C'q{u) and C[ = C[(v) such 
that, for every e G [0, 1] and r > 1, 



C' (r- 2 I {r< i } + ee- r ^I {r> i } ) < F So (Z^ hzts (-r,r) c ) < C[(r~ 2 1 {r< i } + ee~ r ^l [r> i } ) 



Remark. As an immediate consequence of the upper bound of the lemma, we have in dimension d, 
for every r > 0, 

F So (Z (0) hits B(0, r) c ) < C'{ (r + 1)~ 2 (7) 

with a constant C" = C'{(d,iy). 

Proof. It clearly suffices to prove that the stated bounds hold for the quantity ¥$ r (Z^ hits (— oo,0]) 
instead of ¥s (Z^ hits (— r, r) c ). We fix e > 0, and for every x > and t > 0, we set 

q e (x,t) = Fg x (^Z^ does not hit (— oo,0] before time t 

and 

p £ {x) = ¥ Sx (Z {£) hits (-00,0]) = lim f (1 - q(x,t)). 
\ I tfoo 

In this proof only, we write P| for the probability under which £ is a Brownian motion starting from 
x and killed at rate e (upon killing, £ is sent to the cemetery point A and we recall that all functions 
vanish at A). Write 5 := inf{t > : G (— oo,0]}. By standard arguments (see e.g. the proof of 
Proposition II. 3 in |LG99| ). the function q E solves the integral equation 

/•tAS 

q e (x,t)=P £ x (S>t)+E £ x / (r(q E (Z s ,t-s))-q £ (ti s ,t-s))ds 

where we recall that T denotes the generating function of the offspring distribution v. For every 
a 6 [0, 1], set ^(a) = T(l — a) — (1 — a). Note that $(0) = and the function $ is monotone increasing 
under our assumptions. Furthermore, $(a) = \o? + o(a 2 ) when a — > 0. By a monotone passage to 
the limit we get that, for every x > 0, 

p £ (x) + E| [ / *(p e (&)) (fel = P^(5 < oo). (8) 
L Jo J 

It follows that the function p e satisfies the differential equation 

-p" = ep £ + *(p 6 ) 

on (0, oo) with boundary conditions p e (0) = l,p e (oo) = 0. By solving this differential equation, we 
get, for every x > 0, 

du 

- x 



,( X ) yj2eu 2 + 4T(u 

where F(u) = J " $>(v ) dv . Note that there exist positive constants c, d such that cv? < T(u) < c'u 3 
for every u S [0, 1]. The desired bounds then follow from easy analytic arguments. □ 
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2.4 Hitting probabilities for super-Brownian motion 

Let = (Y t {a > 

)t>0 be a super-Brownian motion with branching mechanism i(jr a \ for some a > 0. 
Suppose that starts from under the probability measure P„, for every \i E Mf(R d ). 
The range of Y^ is by definition 

^ (a) )=u fO su pp^ (a) )V 

e >0 \t=e / 

where for every /i E .A/f^B^), supp(/u) denotes the topological support of /z. 

Let D be a domain in M d and let x E D. Consider the process Y^ started from 5 X . We say 
that y (a) hits D c if the range K(Y^) intersects D c . By a famous result of Dynkin |Dy91[ |Dy93| the 
function 

ug ) (x) = -io g (i-p fc (w(yW)nDV0)) . 

is the maximal nonnegative solution of the semilinear partial differential equation ^Au = ipr a -\(u) in 
D. 

Under mild regularity assumptions on D (which hold e.g. when D satisfies an exterior cone condition 
at every point of dD), the function uf\ has boundary value +00 at every point of dD and is the unique 

nonnegative solution of the equation ^Au = ipr a \(u) in D with boundary value +00 everywhere on 
dD. A discussion of this result and related ones can be found in Chapter VI of the book [LG99J. This 
reference considers only the case a = 0, but the same results can be obtained for any a > by similar 
arguments: Note that the Brownian snake approach can be extended from the case a = considered in 
[LG99J to a > 0, simply by replacing the reflecting Brownian motion driving the snake by a reflecting 
Brownian motion with negative drift (see Chapter 4 of I )!.( 102 for a discussion of the snake approach 
to superprocesses with a general branching mechanism). 

We shall be interested in the special case D = A. Recall that E A and that we assume A is a 
domain of class C 2 , meaning that the boundary of A can be represented locally as the graph of a twice 
continuously differentiable function, in a suitable system of coordinates. We write U( a )(z) = u^(x) to 
simplify notation. From the analytic viewpoint, the function ur^ may be constructed as follows. For 
every integer n > 1, let U( a ^ n be the unique nonnegative solution of the nonlinear Dirichlet problem 

f ±Au = ip( a )(u) in A, 
\ u\dA = n . 

Then = lim f U( a ),n as n — )• 00. 

The following lemma records certain analytic properties which will be useful in the forthcoming 
proofs. 

Lemma 2. (i) Let x 6 A. The function a — > u^(x) is continuous and nonincreasing on [0,oo), and 
tends to as a —> 00. 

(ii) For every 5 E (0, dist(0, A c )), let As be the subdomain of A defined as the connected component of 
the open set {x£i: dist(x,A c ) > 5} that contains 0. Then, for every a > 0, ii^(0) tends to u^(0) 
as 5 — > 0. 

Proof, (i) Let us first verify that the function a — > u^{x) is monotone nonincreasing, for every 
x E A. To see this, we apply a standard comparison principle (see e.g. Lemma V.7 in [LG99]) to 
obtain that U( a /) ra < ii( a ) jn if a < a', for every n > 1. It then suffices to let n — > 00. 

Let (ak)k>i be a sequence of nonnegative reals increasing to a E (0, 00). We can set for every x E A 

v{x) = lim I u {ak) {x), 
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and we have v > U/ a -\ . In order to verify that v < Uf a \ , we only need to check that v solves | Av = ipi a \ (v) 
in A (recall that Ur a \ is the maximal nonnegative solution of this equation). To do so, let B be an 
open ball whose closure B is contained in A. For every k > 1, the restriction of ut ak ) to B solves the 
equation \Au = tpu k \(u) i n B. By the probabilistic interpretation of the integral equation associated 
with this PDE (see e.g. Chapter V in [LG99J), this implies that, for every x £ B, 



u (a k ){x) +E a 







TB 



E x[u {ak) (Cr B )] 



where we recall our notation £ for a Brownian motion starting from x under the probability measure 
P x , and tb '■= inf{i > : ^ £ B}. By passing to the limit k — > oo in the previous display, we can 
write 



v(x) + E x 







^(o)(«(Cs))ds 



which is enough to obtain that v solves \Av = ipt a )(v) m B, and therefore in A since S was arbitrary 
Similar arguments show that, if (afc)fc>i is a decreasing sequence of nonnegative reals converging 
to a € [0, oo), then ut ak \(x) converges to Uf a \(x) for every x S A. Finally, the fact that uu-\(x) tends 
to as a —7- oo can be obtained from the comparison principle: If B is a ball such that B C A, the 
restriction of to B is bounded above by the solution vr^ of the linear equation ^Av^ = ao(„) 
in B, with boundary value equal to the restriction of ur^ to B. It is easily seen that Vi a \{x) — > as 
a — > oo, for instance by using the Feynman-Kac formula. 

(ii) Fix a > 0. If < 5 < 8', the closure of A$> is contained in Ag. The restriction of u,\ to yly is 
a nonnegative solution of \Au = "0(a) ( u ) m ^4<5' an d i s thus bounded above by ufy Hence, for every 



(a) 

fixed x S A the function <5 — > u^(x), which is defined for 5 > small enough, is nondecreasing and 



we can set 

v(x) = limu^\(x). 

By the same argument we used to obtain the monotonicity of the mapping 5 — > u^(x), we also have 
v{x) > U( a )(x) for every x S A. To obtain the reverse inequality v < U( a ), it is enough to verify that 
v solves \Av = tfj^(v) in B. But this follows by arguments similar to those we used in the proof of 
part (i) of the lemma. □ 

Lemma 3. For every a > and x 6 A, 

{K(Y^) n A c / 0} = {K(Y {a) ) n A c ^ 0} , P & a.s. 

Proof. The inclusion 

is trivial. To show the reverse inclusion, we may argue as follows. By Theorem IV. 9 in [LG99J (which 
holds under much less stringent assumptions on ^4), the event {TZ(Y^) n A c ^ 0} holds if and only if 
the exit measure of the super-Brownian motion from A is nonzero. Applying the special Markov 
property of superprocesses |Dy93, Theorem 1.3], we see that it is enough to prove that for super- 
Brownian motion starting from a nonzero initial measure supported on dA, the range immediately hits 
(A) c . This is however easy under our regularity assumptions on A. We leave the details to the reader. 
□ 

From now on, we write {Y^ hits F} for the event {TZ(Y^) n F ^ 0}. 
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3 Quenched convergence to super-Brownian motion 

The main goal of this section is to prove Theorem [2] At the end of the section, we also establish 
Proposition [TJ using certain arguments related to the proof of Theorem [2j To simplify notation, we 
set for every w E fi and e E (0, 1), 

v% = VFr ro = |J y/Hxi + Ki). 

iei 

The following lemma identifies a martingale problem solved by our branching Brownian motion 
Z m ' £ . It can be proved by standard arguments (see e.g. Section 9.4 in [EK86J). 

Lemma 4. Let w E fl and e > 0. Under each probability P M , /i E A4 p (M. d ), the process Z ro,£ solves 
the following martingale problem: For every f E (7^(1^) stic/i that < inf / < / < 1, i/ie process 

is a martingale. 

We can derive from Lemma [4] (or from a direct argument) that for every g E C' 2 (R c( ), 

M t (0) := (Zr £ ,5) - (Zo> £ ,9) ~ f (Z7' £ , \Ag - el^g) ds (9) 
is a martingale. An easy computation gives that the square bracket of this martingale is 

[M(g),M(g)] t = f\z^,Vg.Vg) ds + £ (Z^ - Z™f,g) 2 . (10) 

"'° 0<s<t 

The last sum in the right-hand side is an increasing process with compensator 

f{Zr>\(o- 2 + eI T Jg 2 ) ds. (11) 
J o 

The proof of Theorem [2] relies on the following two results, in which we use the notation (3^t)t>o f° r 
the canonical process on Dj^ [0, oo). Recall from (|3|) the definition of the process X m,e in terms 
of Z ro ' £ . 

Lemma 5. For every w E ft and e E (0,1), Zef n ro,£ 6e i/ie Zaio o/ i/ie process X m ' £ under Pr e -ii£ . 
Then, 

(i) For every 5 > and T > 0, f/iere is a compact subset K$,t ofW 1 such that, for every w E fl, 



sup n ro > £ sup y t {K c ST )) < 6 

£6(0,1) V 0<t<T ' y 

(ii) For every g E C l2 (M rf ) andw E fi, i/ie collection of the laws of the process ((34, <?})t>0 under H m,e , 
e E (0, 1), is relatively compact in the space of all probability measures on -Dr[0, oo). 

Consequently, for every w E f2, i/ie collection (n ro ' e ) £g ( ,i) * s relatively compact in the space of all 
probability measures on ( K d)[0,oo). 

The last assertion of the lemma is an immediate consequence of (i) and (ii) using Theorem II. 4.1 
in |Per02| . 
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Proposition 3. Let g G Ci.(R ). There exists a measurable subset ft g of ft such that P(f2 ff ) = 1 and 

the following holds for every w G Q g . For every s,t > 0, for every integer p£N and every choice of 
ti,...,t p e [0,t] and fi,...,f p e C(Mf(R d )), we have 



lim Er„-i 



E-S>0 



Jt 1 > i=i 



0. 



We postpone the proof of Lemma and Proposition [U and explain how Theorem [2] follows from 
these two statements. We choose a countable dense subset G of C\{M. d ) and set 

n' = n *V 

Fix G fi'. By Lemma [5j the collection (n ro ' £ ) eg (o,i) i s relatively compact. Let IT* be a sequential 
limit of this collection as e tends to 0. We deduce from Proposition [3] that, for every g G G, for every 
s, t > and every choice of t±, . . . , t p G [0, t] and f\, ■ ■ ■ , f p G (7 (.M/(M d )) , we have 



IP 



0. 



(12) 



The required passage to the limit under the expectation sign is easily justified: Note that, by comparing 
with the case when there is no killing and using standard results of the theory of critical branching 
processes, we have for every T > 



sup [E [e -i So] sup (X?' £ ,l) 

ee(0,l) V L 0<r<T 



< oo. 



(13) 



Since G is dense in C+(M. d ), another easily justified passage to the limit shows that (1 12 [) holds for every 
g G C+(R ). Thus, II* satisfies the martingale problem for super-Brownian motion with branching 
mechanism i^( K \ as stated in Section |2] Since it is also clear that II* (Yq = 5o) = 1, II* must be the law 
of super-Brownian motion with branching mechanism ift^ started from 5$. This completes the proof 
of Theorem [2] but we still need to establish Lemma [5] and Proposition |3] □ 

Proof of Lemma [5l The compact containment condition (i) in the lemma is immediately obtained 
by observing that X m ' e is dominated by X° = X m '° and by using the convergence of rescaled branching 
Brownian motions (without killing in the obstacles) towards super-Brownian motion. So we only need 
to verify (ii). In the remaining part of the proof, we fix zu G fi and g G Cm.(R ). To simplify notation, 
we set 

X t c = (Xr,g). 

By the remarks following Lemma |4] and an elementary scaling transformation, we have 



Xt = X^ + Mt + V t 



t I 



where 



J" 

Jo 



X^,-Ag-I r ^g) ds 



and M e is a martingale, whose square bracket is given by 

[M e ,M%= f(X^ £ ,eVg.Vg) ds + £ (X?* - X™f,gf 

"'° 0<s<t 
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Furthermore, the oblique bracket of M £ is equal to 

(M £ ,M £ ) t = f (X?> s ,a 2 g 2 + e(Vg.Vg + Irey)) ds. 
Jo 

By standard criteria (see in particular Theorem VI. 4. 13 in Jacod and Shiryaev [ JS87J ) , the tightness 
of the laws of the processes X s , e G (0, 1), will follow if we can verify that the laws of the processes V s 
and (M £ ,M £ ), for e G (0, 1), are C-tight. But this is immediate from the preceding explicit formulas 
and CE). □ 

Proof of Proposition [H Let us fix s, t > 0, t 1} . . . , t p G [0,t] and fi, . . . , f p G C(M f (R d )) . Let 
w G fi and e G (0, 1). Using the facts that T(l) = T'(l) = 1, T"(l) = a 2 and Lemma g] applied to 
the function f(x) = exp { — eg(xy/e)} , we can write 



= E 



'x^ £ ,^{-e 2 Ag + e 3 Vg.Vg) 



+ el r ^(l-e- £9 )e e9 + e e9 (T(e- E9 ) - e~ £9 )\ e'^^du i f[ fi(X™' £ ) 



1=1 



i=i 



(14) 



where ||r/ e || < C\(g)s for some constant C\(g) depending only on g, and ||ry e || — >• as e — >■ 0. On the 
one hand, using Fubini's theorem, the fact that {{X™ ,£ , l)) r >o is a supermartingale and the inequality 
Efe-iiaJCXj 7 ' 6 ,!)] < 1, we have 



E 



ti J 



< (1(1^^11+2(7! (5))+Yll%l|)(nil^ll) E [- 1 15o[/ +S W' £ ) l 
<(C 2 ( 5 ) £ + (a 2 /2)||%||) f ni|/i||^ 



i=l 



with a constant C2{g) depending only on g. On the other hand, 



E [e- 1 ]5o 



Jt J i=l 

t+ Vr e ,^> e-^'.ddu} n/ ; (in 

j i=i 

1/4 ^ i=l 



The absolute value of the second term in the right-hand side of ([15]) is bounded above by 



(15) 



i=i 



E [e- 1 ]5o 



(X^I^ -«)<?) | 



(in. 
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Consequently, going back to (|14p and using the preceding estimates, we obtain for every w G fi and 
£€(0,1), 



E[.-i]«b 



< 



n 

i=l 



{C 2 ( 5 )£ + (a 2 /2)||77 e ||}s + r £ (tu, 5 ,t,t + s) 



(16) 



where r e (w, g, t\, i 2 ) is defined for < t± < i 2 by 
Lemma 6. Let u > 0. TTien, P-a.s. 



(w, g, h,t 2 )= ^ E [£ _ 1]5o [|<A^, (ir* - k) 5 ) | 
Ju L 



<iu. 



(17) 



Assuming that the lemma is proved, we can readily complete the proof of Proposition [3] Using 
Fubini's theorem, we may find a set fl g with P[fi g ] = 1, such that, for every w G fi 9 , the convergence 
in (|17p holds simultaneously for all u > 0, except possibly on a set of zero Lebesgue measure (depending 
on w). Since for every w G f2, e 6 (0, 1) and u > 



E 



(AT' £ , (Ir«,-K)5>| <|H|E [£ - 1](5o [(X-'M)]<|| 5 || 



dominated convergence guarantees that for each w G f2 9 , lim e _-.o r e( ro > 9^^ + s) =0 for all £, s > 0. It 
follows that the right-hand side of (jl 6[) tends to as e — > when G fi 9 , which completes the proof 
of Proposition [3] □ 

Proof of Lemma [6l Recall that £ ro ' e denotes standard d-dimensional Brownian motion killed at 
rate e within Y m . We also use the notation x™' 6 101 Brownian motion killed at rate 1 within F^. Both 
processes and x ro ' e start from x under the probability measure P^. 

We first recall classical moment formulas for branching Brownian motion. For every x G M. d , k G N, 
and every bounded measurable function h on M. d , we have 



E k5x [{Zr,h)]=kE x [hti7' £ )) 



and 



E kSx [(zr,h) 2 } = kE x [h(^' E f] + Kk - i)E x [h(^nr + ^e x 



dsE if 4h(^l)f 



These formulas are easily derived, first for k = 1, from the well-known formula for the Laplace functional 
of Z™' £ under Wg x : See e.g. Proposition II. 3 in [LG99J. 

Recalling the definition of X m,e in terms of we get similar formulas for the first and second 

moment of {Xf' £ , h). In particular, for s > and for every w G fi and e £ (0,1), we have 



E 



<AT' £ , (l a - K ) 5 > 2 1 = e 2 ^ 1 ] 2 - [e-^A^^s^gf +e[e- 1 ]A^{Q,s,g), (18) 



where, for every x G 



£ (z, S , 3 ) := E x [l^Jxf' e ) " «)s(xr ,e ) 
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and 



Af' e (x,s,g) := a 2 E a 



v,g) dv 



We claim that, for every x G and s > 0, 



P— a.s., 



lim E a . 



(ii*(x^-«Mxr) 



0. 



(19) 



Assume for the moment that the claim holds. By taking x = in (I19p . we obtain that A™' £ (0, s, g) 
tends to as e — > 0, P-a.s. Consider next A™' 6 (0, s, g). The second term in the formula for A™' e (x, s, g) 
obviously tends to uniformly in x, s and independently of the obstacles, as e — > 0. To handle the 
first term, use Fubini's theorem to obtain that, for every w belonging to a set ft of full probability, 
there exists a set M w C M. d x M + of zero Lebesgue measure such that the convergence in (1 19[) holds 
simultaneously for all (x, s) ^ M m . We can then write the first term in the formula for A™' £ (0, s,g) as 
follows: 

2 



dv 



o 



dyp^(0,y) E y (H(xr4)-*)S(X^) 



(20) 



where p™' s (-, •) is the transition kernel of x m ' e a t time v. Plainly, we have p™' £ (x, y) < p v (x, y) for each 
v > and e G (0,1), where p v (-,-) is the transition kernel of standard Brownian motion at time v. 
Using this bound and dominated convergence gives us that the quantity in (120p tends to as e — > 0, 
for every w G ft. Hence A™' (0, s, g) tends to as e — >• 0, P-a.s. The result of Lemma [6] now follows 
from (|18p and the Cauchy-Schwarz inequality. 

We still have to prove our claim (jl9j) . We fix x G R d and s > 0. Let G f2 and e G (0, 1). From 



the definition of Brownian motion killed at rate 1 in T^, we have 



E, 



s ) 



k)q{xT £ ) 



E, 



(Irj, (6) " «) <7(6) exp { - J' lr^u)dv] 



(21) 



Let r] > and choose G (0, s) such that e 9 > 1 — r\. By the Markov property applied at time s — 6, 
the quantities in (I2ip are equal to 



E, 



exp { - J H (£«)<*«} E 5s „ e [(l a (&) - «)(/(&) exp { - jf 
6XP { ~ ii Ir|> ^"^ n } 



E x 
+E 



exp | - ^ Ir*, E s 



ir^) 

(lr*,(&) exp 



(22) 



i 



The condition 1 — e e < T] entails that the absolute value of the second term in the right-hand side of 
(122p is bounded above by \\g\\rj. Suppose we know that for every y G M rf we have 



P— a.s., 



lim E 

e^0 



(lrg,(^)-«)^) 



0. 



(23) 



Then, using the fact that the law of i s -o is absolutely continuous with respect to Lebesgue measure 
(and Fubini's theorem to see that the convergence in (I23p holds simultaneously for almost all y G M. d , 
P-a.s.), we conclude that the first term in the right-hand side of 022 p converges to with P-probability 
1 as e — > 0. Consequently, we have P-a.s. 



lim sup 

e^0 



E x 



(Ircjx7' e )-K)g(xl 



Since r\ was arbitrary, our claim (I19p follows. 
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It thus remains to prove that (|23p holds. We fix y £ M d and 9 > 0. In the following, £' stands for 
another Brownian motion independent of £, which also starts from y under the probability measure 
P y . For each e £ (0, 1), we have 



E 



Ej(l r57 (&)-K)s(&)] : 



E 

E, 



(24) 



where the last line uses Fubini's theorem and the definition of k. Recall that the measure is 
supported on compact sets which are contained in the fixed ball B(0,ro). If \x — x'\ > 2ro, the sets 
{(z, K) : x e z + K and K C B(0, r )} and {(z, K) : x' £ z + K and K C B(0, r )} are disjoint, and 
so the events {x G and {x 1 £ are independent under P. Recalling that = y/e T w , we see 
that if \x — x'\ > 2ro"v/e we have P[x £ T^; x' £ T^] = P[x £ P^jPfa:' G = k 2 , which enables us 
to write 



Going back to (|24p . we obtain 



E 



H*[fe(6)-«)s(60] : 



E,, 



= IM| 2 Po[|£ 2 e| <2r ov ^] 



« 2 } 



0d/2 



T d/2 



(25) 



where the constant C > depends only on rg. Let t/ > 0. By the Markov inequality, we can write 



|E,[(Irs,(6) -*)»({»)] I >v 



if 



< C- 



2Qd/2 



T d/2 



Applying the last bound with e = £/% = k 3 , for every k £ N, yields a convergent series, and so by the 
Borel-Cantelli lemma we have P-a.s., 



lim sup 

k— >oo 



< 



Since rj > was arbitrary, the convergence (|23p holds along the sequence (£k)k>i- 
To complete the proof, we set for every w £ ft and e £ (0, 1) 

EW = E„[ll- (&)!/(&)] 
(recall that y and are fixed). We shall prove the following result: For every w € ft, 



lim sup 

fc ^°° (fc+l)-»<e<fc-3 



17, 



0. 



(26) 



Combining (|26p with the fact that the convergence (|23p is true along the sequence (sk)k>i w iU then 
lead to the desired conclusion. In order to prove (|26p . we first note that for every e £ (0, 1), 



l/tz 



^r%(x) g(x)e 2a 2 dx 



Ir s (x) h(x)dx, 



(27T0) d / 2 J Rd " 

with a function /i £ C?.(R ) which depends on y and 0, but not on tu. Furthermore, for any fixed 
7] > we can find a large closed ball B centered at the origin and such that J Bc h{x)dx < rj. Hence, if 



wc set 



u' 



B 



h(x) lye (x)dx, 
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we have |C/ roj£ — U' w e | < r\ for every e £ (0, 1). Thanks to this remark, it is enough to prove that (|26p 



holds when U miE and U^ tEk are replaced by U!^ e and U' VJfih respectively. 
Let k G N and e € [(& + 1)~ 3 , A; -3 ). We have 



■ e 



d/2 
d/2 



•1/2B 



h[x\/e)lr m (x)dx 



k 3 / 2 B 



d/2 



{£~ 1 / 2 B)\(k 3 / 2 B) 



/i(x-v/e)lr ra (^)^- 



(27) 



Now, the first term in the right-hand side of (|27p is equal to 



d/2 1 



k 3d/2 



k 3/2 B 



k 3/2 



l T ^{x)dx + e d/2 



k 3/2 B 



h(x\/e) — h 



k 3 / 2 



Ir^(x)dx 



(sk 3 ) d/2 U^ £k + i(xu,s,k), 



where 



i(zu,e,k) = (ek 3 ) d / 2 / \h{{k^ 2 ^-e)y) - h{y)] I r ^(k 3 / 2 y)dy. 



B 



Note that < 1 — k^^yfe < with a constant C independent of e and k, from which it easily follows 



that 



sup 



a 



\l(w,e, k)\ < — ||V/i||, 

(fc+l)- 3 <e<fc- 3 k 

with a constant C depending only on B. 

Similarly, the second term in the right-hand side of (I27p is bounded above by 



e d ' 2 \\h\\[e- d l 2 -k^ 2 ]\ d {B) 
Finally, from (l27|l and the preceding estimates, we have 



< 



C" 



U'\<(l-(sk 3 ) d / 2 )U' + 



and the convergence f)26[) follows. This completes the proof. 



C'\\Vh\\ C" 



k 



+ 



A: 



□ 



Proof of Proposition [TJ Clearly we get a lower bound for the quantity Fs (Z w ' £ hits B(0, R) c ) if we 
replace by branching Brownian motion killed homogeneously over M d at rate e. Part (i) of the 

proposition thus follows from the lower bound in Lemma [1] 

Let us turn to the proof of the first assertion in (ii). We have the following inequality: For every 
zu eft and e G (0, 1), 



5o [Z™> £ hits B(0, R) c ] < P [C ,£ hits B(0, R) c ] . 



(28) 



Indeed, using the formalism of Subsection 12. 3| the criticality of the offspring distribution can be used 
to check that the right-hand side of (128[) is just the expected value of the number of those historical 
paths u v that first exit B(0,R) during the interval [a v ,(3 v A C?' £ ). Alternatively, it is easy to derive 
an integral equation similar to (j8j) for the function x — > Fg x [Z ro ' £ hits -B(0, R) c ] , and the bound ( 128)) 
then trivially follows from this integral equation. 

Then, let us bound P ® P [£ ro,£ hits B(0, R) c ] . We write £ ro ' £ ' fc for the fc-th coordinate of for 
every k = 1, . . . ,d. First, observe that 



P <g> P [C ro,£ hits B(0, R) c ] < P ® Po [C ro,£ ' fc hits (-R/Vd, R/Vd) c ] . 



(29) 



k=l 
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Clearly, we can restrict our attention to the first term in the sum. Define N E r := [Ry/e/y/d\ and for 
every j G Z, set T? = mf{t > : = J£~ 1//2 }, where £ x stands for the first coordinate of £. In this 
notation, we can write 

P® Pof^' 1 hits (-R/Vd, R/Vd) c ] 

r' 6 ' 1 hits(-iV £ ^- 1/2 ,iV £ ^- 1/2 ) 



<P®P 
<E®E 



exp < — e 



+ E®E 



exp<^ -e / e ' Ir w (6)rfs 



(30) 



Consider the first term in the right-hand side of (j30[) . Recall the definition ([I]) of the set of obstacles, 
and set for every j G N and e > 0, 



I':={iel: Xi G ((j - l)e' 1/2 ,je 



1/2 ,---1/2^ 



»d-l 



r ro (j,e) := |J (xi + Ki) and := y/i T m (j,e). 



Note that the random sets r ro (j, e), j G N, are independent under P, by properties of Poisson measures. 
We have then 



E<g>E 







exp | 









Ir w (C.)d* 



<E® E 



^EtgiEn 



N, 



exp < — e 



I 



exp<^ -e / 1^(1,6) (6 



(31) 



where the equality comes from an application of the strong Markov property of £, together with the 
independence of the random sets T ro (j, e) and the fact that the distribution of each of these random 
sets is invariant under translations by elements of {0} x By scaling, if T\ denotes the entrance 

time of £ into [1, oo) x we can write 



a £ := E® E 
We then observe that 



exp | 


-/ 




Jo 



E(g)En 



exp 



Ti 



Irg,(i)(6)ds 



Ti 



Hd) (6)1(0,1)^)^ « / Tl i (0) i)(e s V 

£ ^ u JO 



(32) 



where the notation — ?■ refers to convergence in probability under P(g)Po. To see this, we use arguments 

similar to the proof of Lemma El Notice that P[y G 1^(1)] < k if y G (0, 1) x with equality if 

y G (ro^/e, 1 — roy/e) x Using this remark, and the same argument as in the derivation of (125 p . 

we can write for every fixed u > 0, 



E® Ei 



TiAu 







E<g>E 

: Eq 



JO 
TiAm /-TiAu 



(%(i)(6) - «)(% ( i)(6) - k) I ( o,i 3 (£) Vi) $)dsdt 



— Eo 
<E 



ii Jo 

TiAu /-TiAu 



{p[6,& g r e ro (i)] - « 2 } 1(0,1)^) Vd^* 1 )^ 



+ 0(e 



l/2\ 



JO 

u fU 



h\t'-tt\<2r ,M ho,i) \o,i) dsdt 



+ 0(e 



1/2 > 



JO 



+ 0(e 1 /2 ) 
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where the error term 0(e 1 ^ 2 ) corresponds to the contribution of times s,t such that or belongs 
to the set (0,roy/e\ U [1 — ro^/e, 1). The preceding quantity tends to as e — > 0, which yields the 
convergence (j32~j) . Since the limiting variable in (l32~j) is (strictly) positive a.s., we can find e\ G (0, 1) 
and ci < 1 such that a e < c\ for every £6 (0,£i). Using (|29p . (|30p and (|3ip . we arrive at 

P®P [e ro ' £ hits B(0,R) c ] < 2dc^ e ' R , 

for e G (0, ei). This completes the proof of the first assertion in (ii). 

Let us turn to the second assertion. Let C2 > be a positive constant whose choice will be specified 
later. By simple comparison arguments, it is enough to prove the desired estimate when R is of the 
form R = 2 k , for k G N large enough, and e is of the form e = 2 _J , with j G {0, 1, . . . , 2k} such that 
e > C 2 (loglog R) 2 /R 2 . 

By the first assertion in (ii) and the Markov inequality, 

< exp(-Cii? v / e/2). (33) 
However, if R = 2 k and e > C 2 {\og\ogR) 2 / R 2 , we have 

Ry/e > VC~2 log log R = ^(log k + log log 2). 

Using this bound, we can choose the constant C2 sufficiently large so that we get a convergence series 
when we sum the right-hand side of H33[) over all R = 2 k and e = 2~- J for j G {0, 1, . . . , 2k} such that 
e > C 2 (loglog J R) 2 /i? 2 . The Borel-Cantelli lemma now yields the desired result. 

□ 



P 5o [Z ro ' e hits B(0,R) c ] > exp(-Cii? v / e/2) 



4 Proof of the main result 

In this section, we prove Theorem [TJ We fix the environment w such that the weak convergence of 
Theorem [2] holds, and derive the convergence in Theorem [1] for this fixed value of the environment. For 
the sake of simplicity, we shall omit w in the notation and write Z e instead of Z ro,£ , and X s instead 
of X™' £ . 

We shall verify that for any increasing sequence {R n )n>i of positive reals converging to +00 and 
any sequence (e n )„>i of nonnegative reals such that £ n Rn [0, 00], we have 

lim R 2 n F So {Z £ - hits R n A c ) = ur Ka) (0), (34) 

where U(oo)(0) = by convention. 

The statement of Theorem [1] follows from this convergence. Indeed, if the conclusion of the theorem 
fails, then we can find a sequence R n f 00 and a sequence (e n ) of nonnegative reals such that, for every 
n > 1, 

\R 2 n Ps Q (Z £n hits R n A c ) - u {KenK) (0)\ > 6 

for some constant 5 > 0. By extracting a subsequence, we may assume that £ n R\ — ► a, G [0, 00] and 
thus obtain a contradiction with (j34)) since we know from Lemma [2] that the mapping b — > un,\ (0) is 
continuous on [0,oo]. 

In proving (|34p . we may assume that e n — > as n — > 00. Indeed, suppose that 034 p holds in this 
particular case and let (e' n ) be a sequence that does not converge to 0. If the sequence e^-R 2 converges 
then necessarily its limit is +00, and we can find another sequence such that < e 1 ^ < e' n , e" — > 
and e'nRn — > 00. So, if we know that (j34j) holds in the case when the sequence (e n ) tends to 0, we 
obtain 

lim RlF So {Z 6 "* hits R n A c ) = 0. 
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However, from the inequality e" n < e' n and a coupling argument (obvious if one uses the construction 
described in Subsection I2.3p . we get the same result for the sequence (s' n ). 

A similar comparison argument shows that it is enough to prove f|34[) in the case when a < oo. 
Otherwise, it suffices to replace e n by e n A bR~ 2 and let b — > oo, using the fact that un,)(fi) — )• as 
b — >• oo. 

Let us now proceed to the proof of (|34p . We fix the sequences R n f oo and e n — > such that 
e n R\ — > a £ [0,oo). We first assume that a > 0. The case a = will be discussed at the end of the 
section. 

Let B be a closed subset of R d . For every e > 0, we have by the definition of X s 

P [£ - 1]5o (X £ hits B)=¥ [e - 1]5o (3t > : Xf(B) > 0) 

= F [e -i ]So fat > : J Z^ H (dx)I B (x^I) > 0\ 

= P [E - 1]5o (Z £ hits e-^B) 

= 1 - F 5o (Z £ does not hit e~ 1/2 B) [e ~ 1] , 

since (for a fixed environment) the law of Z e under P[ e -i],5 is obtained by adding [e _1 ] independent 

1/2 

copies of Z £ under P,5 . Applying the preceding identity with e = e n and B = b n A c , where b n = e n R n , 
gives us that 

1 - P SQ (Z £n does not hit R n A c )^ = F [£ -i ]So (X 6n hits b n A c ) (35) 

By Theorem [21 we know that the law of X 6n under P^-ij^ converges as n — > 00 to the law of super- 
Brownian motion with branching mechanism ipr K \ started from ^o- The next lemma is essentially a 
consequence of this convergence. We use the notation of Subsection 12.41 

Lemma 7. We have 

lim F,-i ]s (X £ " hits b n A c ) = P So (yW hits bA% 

where b = yfa = lim b n . 

We postpone the proof of Lemma [7] and proceed to the proof of (|34p . in the case when a > 0. By 
the results recalled in Subsection 12.41 we know that 



P 5o (y (K) hits bA c ) = 1 - exp(-u(0)) , (36) 

where the function (v(x),x S bA) is the unique nonnegative solution of the singular boundary value 
problem 

J ±Av = ip( K )(u) in bA, 
\ U\d(bA) = +OO • 

It is immediate to verify that U( Ka )(x) = av(bx) for every x G A, and in particular U( Ka )(0) = av(0). 
From 035 p . (|36p and Lemma [71 we obtain 

lim (1 -F 5o (Z 6n hits R n A c )) [e " 1] = exp(-u(0)) 

n— >oo 

and thus 

lim e^FsJZ 6 " hits R n A c ) = v(0), 

n— >-oo 

or equivalently, since e n R\ — > a, 

lim R 2 n ¥ So {Z £n hits R n A c ) = av{0) = u, Ka) (0). 
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This completes the proof of (j34[) . in the case a > 0. □ 
Proof of Lemma By replacing A with bA, we may and shall assume in this proof that 6=1. We 
thus have b n — > 1 as n — > oo. We first prove that 

liminf P r -i u (X £n hits b n A c ) > P S JY M hits A c ). (37) 

n— >oo l £n J 00 

By Lemma |3l the events {yw hits >l c } and {Y^ hits ^4 C } coincide a.s. We can then find a countable 
collection (</?i)j>i of continuous functions with compact support contained in (A)°, such that 

{KM hits (i) c } = { sup (sup (Y t {K \^)) > o), P tfo a.s.. 

i>l V t>0 7 J 

Hence, if {tj)j>\ is a sequence dense in [0,oo), we have 

P So (Y^ hits (A) c ) = Jim t Ps ( sup ( sup (Y^,^)) > o) . (38) 

N^oo \ i<j<iv v l<j<N J J 



However, Theorem [2] implies that, for every N > 1, 



v l<i<iV v l<j<N v l<i<W v l<j<N 



liminf P r sup ( sup (Xf»,^)) > 0) > P S J sup ( sup (Y t M , c^> ) > ) . (39) 



Recall that b n — > 1, and note that the support of each function fi is at a strictly positive distance of 
the set A. As a consequence, for every fixed N, the support of ifi will be contained in b n A° for every 
£ = 1, . . . , N, as soon as n is large enough. Hence, for all large enough n, 

( sup ( sup (X% n ,<pi)) > o) C {X £n hits b n A c }. 
^ l<i<JV ^ l<j<N 3 ' > 

Using this inclusion and then (|39p and (|38p . we immediately obtain (j37J). 
We next turn to the more difficult upper bound 

limsup Pr -i,. {X £n hits b n A c ) < P 5o {Y {K) hits A c ). (40) 

We fix 5 > small enough so that the closed ball of radius 45 centered at is contained in A. As in 
Lemma [2j we let A$ be the connected component of the open set 

{x G A : dist(x,^ c ) > 5} 

that contains 0. We denote the exit measure from As for the rescaled branching Brownian motion X Sn 
by £ ?. In other words, the measure £^ is equal to e n times the sum of the Dirac point masses at all 
points of OAs which are first exit points from As for one of the historical paths associated with X £n 
(these historical paths are defined in Subsection 12.31 for the branching Brownian motion Z £n , and this 
definition is extended to X £n by an obvious scaling transformation). 

Let $ be a continuous function on R rf such that 0<<3?<l,<l ) = 0on A^g and $ = 1 on A^fi. Then, 
for every r] > and p > 0, 



P [E -i ](5o (A-" hits b n A c ) =F [£ - 1]So (X £ " hits b n A c , (SI 1 , 1) < r/) 

+ F [£ - 1]S (X £ " hits b n A c , > V, J o (Xl n ,®)ds < p 

+ P [£ - 1]5o (V" hits b n A c , (£f,l) > r,, J™(X £ s -,<S>)ds > p\ (41) 

Let a n (f]), /3 n (rj, p) and "fniViP) be the three terms appearing in the right-hand side of (|4ip in this 
order. 
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We first bound a n (jj). Provided n is sufficiently large, b n A c is contained in Ag, 2 and thus 

a n ( V ) < ^ ]So ( X£n hits A 6/2> ^ < 

Note that the times at which the historical paths of X £n exit Ag form a stopping line in the sense of 
[Ch9lj. We can thus apply the strong Markov property at a stopping line (Proposition 2.1 in [Ch9fj) to 
see that a n (5) is bounded above by the probability for a branching Brownian motion (without killing) 
starting initially with less than rje^ particles, that one of the historical paths reaches a distance greater 
than 5/(2yfe^) from its starting point (to be precise we need a slight extension of the results in |Ch91| . 
since our spatial motion is not standard Brownian motion, but Brownian motion killed inside 
The estimate ([7J now gives 

a n (ri)<C'{(d,v)^. (42) 



Then, we have 



Recall that $ = 1 on A%g and in particular $ = 1 on B{x,5) for every x € dAg. We use the strong 
Markov property at the same stopping line as in the previous argument, together with a simple coupling 
argument, to write that 



where X £n is defined in terms of a branching Brownian motion Z £n in the same way as X £n was 
defined from Z £n . This branching Brownian motion Z £n has the same offspring distribution as Z £n , 
but particles are now killed at rate e n homogeneously over M. d . Furthermore, Z £n also starts from Mo 
under the probability measure Pfc<5 - 

By Proposition [2j the law of (X^ n )t>o under Pr^-ii^ converges as n — > oo to the law of KW under 

P^Sq (in the notation of Subsection I2.4[) . Noting that, for every fixed s > 0, Yg a.s. does not charge 
the boundary of the ball B(0, 5), it follows that 

The continuity of sample paths of Y^ 1 ' ensures that Poa(r], p) — > as p — > 0, for every fixed r\ > 0. 
For the term "f n (V: p)> we simply use the bound 



oo 



7n(r?, p) < F [e -i ]So I J (X £ s -, <£}ds > p 
This bound and the weak convergence of Theorem [2] imply that 

limsup 7 n(??,p) < Pj f°° {Yj K \$)d8 >p)< P So (yW hits A 



since $ = on A$g. (To justify the first inequality in the last display, we also use the fact that the 
extinction times of X £n under P^-ij^ are stochastically bounded, which follows from a standard result 
in the case without killing.) 

To complete the argument, fix i? > 0. By Lemma [2] (ii), we can choose 5 > sufficiently small so 
that 

P So (Y^ hits A% & ) < P So (yW hits A c ) + f . 
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From 042 p . we can then choose r] > sufficiently small so that for all large n, 

OLnWl < -. 

Finally we choose p > such that /3 o('?)P) < f ■ From (|4ip and the previous estimates, we obtain 
limsup P fe -i la {X £n hits b n A c ) < P So {Y^ hits A c ) + 0, 

and since i? was arbitrary this completes the proof of (|40p and Lemma [7] □ 
We still have to discuss the case a = in (1 34 p . So, let us consider two sequences (e n )n>i and 
(Rn)n>i such that e n -^n ~~ ^ 0- Let ap > and = e ra V (aoi?~ 2 ). Since e n < e' n , we have 



liminf Rlfs n {Z £n hits i? n yl c ) > lim inf 22? F Sn (Z £ * hits R n A c ) = U( Kan) (0), 

n— »oo n— >oo 



by the case a > 0. By Lemma |2](i), ii( Kao )(0) can be made arbitrarily close to tt( O )(0) when a$ is small, 
and so 

liminf R 2 n W 5 JZ £n hits R n A c ) > u (0) (0). 



To obtain the corresponding upper bound, a similar coupling argument shows that it suffices to consider 
the case when e n = for every n, that is when there is no killing inside the obstacles. Hence, consider 
the branching Brownian motion Z° = Z m, ° (the notation is even more legitimate since does not 

depend on zu). For every p > 0, define a rescaled version of Z° by setting 

(X i f\v)=pf Z%(dx)v{p 1/2 x). 

By Proposition [2j the law of (x[ p ^)t>o under Pr^-ii^ converges to the law of under Pg as p tends 
to 0. Set p n = 22~ 2 , in such a way that 

{Z° hits R n A c } = {X {pn) hits A c }. (43) 

A simplified version of the arguments of the proof of Lemma [7J shows that 

lim sup Pr -i ]s (X ipn) hits A c ) < P 5o {Y (0) hits A c ) = 1 - exp(-u (0) (0)). 

n— >oo 

Arguing as in the first part of the proof of the theorem and using (J43j) yields 

limsup^P 5o (^ £ ™ hits R n A c ) < u (0) (0), 

n— >oo 

which completes the proof of Theorem [T] □ 
Proof of Corollary [H For every r > 0, 



W\/^ W,e < r) = (1 - P 5o (Z ro - £ hits 5(0, e-^r)^ 



It- 



However, Theorem [T] shows that, ~P{dm) a.s., 



2 

n e P 5o (Z ro ' £ hits B^e-VV) ) = ^ x (^P 5o (^ ro ' e hits 5(0, £- 1/2 r) c ; 
converges to w° r 2(0) as e — > 0. The desired result follows. □ 
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Abstract 

We derive asymptotics for the quenched probability that a critical branching Brownian motion 
killed at a small rate e in Poissonian obstacles exits a large domain. Results are formulated in terms 
of the solution to a semilinear partial differential equation with singular boundary conditions. The 
proofs depend on a quenched homogenization theorem for branching Brownian motion among soft 
obstacles. 

Keywords and phrases. Branching Brownian motion, Poissonian obstacles, super-Brownian 
motion, escape probability, homogenization, semilinear partial differential equation. 
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1 Introduction 



In the present work, we are interested in the long-term behaviour of branching Brownian motion killed 
in Poissonian obstacles. Let us start by describing a simple special case of our results. We consider 
a critical branching Brownian motion in M, d (d > 1), where all initial particles start from the origin. 
£T) ■ We assume that particles are killed at a (small) rate e > within random balls of fixed radius, whose 
, centers are distributed according to a homogeneous Poisson point process on M. d . Then, how many 
initial particles do we need so that, with high probability, one of their descendants reaches distance R 
from the origin ? Let p £ (R) be the (quenched) probability for our randomly killed branching Brownian 
motion starting with a single particle at to visit the complement of a large ball of radius R centered 
at the origin. The preceding question is equivalent to determining the limiting behaviour of p £ (R) 
^ ■ when e tends to and simultaneously R tends to infinity. 

The answer involves several regimes depending on the respective values of e and R. If e is small 
in comparison with 1/R 2 , the killing phenomenon does not matter and the result is the same as if 
there were no killing: p e (R) behaves like a constant times 1/R 2 (informally, the branching process 
must survive up to a time of order R 2 so that at least one of the particles travels a distance R, and 
well-known estimates for critical branching processes then lead to the correct asymptotics). On the 
other hand, if e is large in comparison with 1/R 2 , then the probability p £ {R) decreases exponentially 
fast as a function of Ry/e: See Proposition [1] below. 

Our main results focus on the critical regime where eR 2 converges to a constant a > 0. We show that 
the probability p £ {R) behaves like R~ 2 , as in the case without killing, but with a multiplicative constant 
which depends on a and can be identified as the value at the origin of the solution of a semilinear partial 
differential equation with singular boundary conditions. A key tool to derive these asymptotics is a 
quenched homogenization theorem which shows that our branching Brownian motions among obstacles, 
suitably rescaled, are close to super-Brownian motion killed at a certain rate depending on a. 

Let us formulate our assumptions more precisely in order to state our results. First, let us define 
the collection of obstacles. We denote the set of all compact subsets of M. d by K. This set is equipped 
with the usual Hausdorff metric djj- Recall that (/C, dn) is a Polish space. For every r > 0, fC r denotes 
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the subset of K, which consists of all compact sets that are contained in the closed ball of radius r 
centered at the origin. Let G be a finite measure on /C, and assume that G is supported on /C ro for 
some r > 0. Let 

iei 

be a Poisson point measure on ~K d x K, with intensity A^ <g> G, where A^ stands for Lebesgue measure 
on M d . We assume that this point measure is defined on a probability space (ft, P) and we denote the 
generic element of fl by w. Our set of obstacles is then defined by 



T zu = \J(x l + K i ), (1) 



iei 

where obviously X{ + K-i = {z = Xi + y : y G Ki}. Note that we use the notation T ro to emphasize 
that the set of obstacles depends on the variable w representing the environment. Let us also define a 
constant n by 

k = P(0 G r w ) = 1 - exp ( - y G(dif) Ad(A-)) . 

To avoid trivial cases, we assume that re > 0, or equivalently G(Ad(-fQ > 0) > 0. By translation 
invariance, we also have P(x G = re for every x 6 R d . 

Let us now introduce the sequence of branching Brownian motions of interest. Given zu £ £1 and 
a parameter e > 0, we consider a branching Brownian motion on R d such that 

• each particle moves around in W 1 according to the law of Brownian motion killed at rate e within 
r • 

L ZD 1 

• each particle branches at rate 1. During a branching event, the particle generates a random 
number of offspring, according to an offspring distribution u which has mean one and finite 
variance a 2 £ (0, oo). 

This branching Brownian motion is denoted by Z ro ' £ = (Z^ u ' £ )t>o, where Z^' e stands for the sum of 
the Dirac point masses at the particles alive at time t. The processes Z™' £ are defined on a probability 
space Q. For every finite point measure fi on M. d , we use the notation P M for the probability measure 
on under which each of the processes Z w ' £ starts from fi. 

Let A be a bounded domain of class C 2 in M. d containing 0. We say that the branching Brownian 
motion Z™' £ hits A c if there exists t > such that Z™' £ (A C ) > 0. We are interested in asymptotics for 
the quantity 

F 5o (Z m ' £ hits RA C ) 

when R — > oo and e — > 0. Here we use the obvious notation RA = {z = Ry : y £ A}. 

Theorem 1. For every a > 0, let U( a ) = (u^(x),x G A) be the unique nonnegative solution of the 
singular boundary value problem 

1 a cj2 2 

-An = — u +au in A, ... 

2 2 (2) 

, u\dA = +oo • 
Then, 

lim ( sup R 2 F So {Z m ' £ hits RA C ) - u (Kem (0) ) = , P(dw) a.s. 

R-^oo \ £>0 J 
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Remark. The existence and uniqueness of the nonnegative solution of the problem (|2|) is a special 
case of Theorem 4.12 in |Ver96j . Alternatively, a probabilistic proof can be given along the lines of 
[LG991 Chapter VI] 

Let us state a corollary of the theorem, which is motivated by the simple question we asked at 
the beginning of this introduction. For every a > 0, we denote by u°, the solution of the boundary 

problem (|2]) when A is the open unit ball of IR d . 

Corollary 1. For every e £ (0,1), let n e be a positive integer. Assume that en e — > b as e — > 0, 
for some b > 0. Denote by R m ' £ the maximal distance from the origin attained by a particle of the 
branching Brownian motion Z w,£ . Then, P-a.s., the law of yfe R m ' £ under Pn E <5 converges as e — > 
towards the probability measure ir^ on M + defined by 

7r 6 ([0,r]) = exp | - ^^ Kr2) (0)| 

for every r > 0. 

In the setting of Theorem [T] it is not hard to see that ui a \ (0) tends to as a — > oo (see Lemma 
[2] below) and thus this theorem does not give much information when R — > oo and e — > in such a 
way that eR 2 tends to oo. In that case, the next proposition provides an exponential decay, which 
contrasts with the preceding theorem. Since our bounds are clearly not optimal, we consider only the 
case when A is a ball. We denote the open ball of radius r centered at the origin by 5(0, r). In the 
general case we may apply the bounds (i) and (ii) of the proposition after replacing A by a ball B(0, r) 
such that -B(0, r) D A or B(0,r) C A respectively. 

Proposition 1. Recall that v denotes the offspring distribution of the particles. 

(i) There exists a positive constant Co = Cq{v) such that, for every w G ft, R > 1 and e £ [1/R 2 , 1], 

¥ 5o (Z m ' £ hits B(0,R) c ) > C e exp(- J R V / 2e). 

(ii) There exists a positive constant C\ = C\(d, v, 0) such that for every R > 1 and e G [0, 1], 

P ® F 5o (Z™' £ hits B(0, R) c ) < exp(-CiiVe)- 

Moreover, we can find two positive constants C% = Ci(d, v, O) and C3 = C^{d, is, O) such that P-a.s., 
for every sufficiently large R and every e G [C2 (log log R) 2 /R 2 , 1], 

F 5o (Z m ' £ hits B(0,R) c ) < exp(-C 3 iVe). 

Part (i) of the proposition is derived from an estimate about branching Brownian motion killed 
homogeneously at rate e, which explains why this bound holds for every w £ fi and does not depend 
on the measure 0. The bounds in (ii) follow from an estimate for Brownian motion killed in soft 
obstacles. The first assertion in (ii) may be compared to Proposition 5.2.8 in [Szn98j. 

Proposition [T] only gives rather crude estimates, and it would be of interest to obtain more precise 
information on the decay of the quenched probabilities Fg (Z m ' £ hits RA C ) in the case when eR 2 tends 
to 00. Because our branching mechanism is critical, the following upper bound holds for any w £ ft, 
e > and R > : 

P 5o (Z m ' £ hits B(0, R) c ) < E 5o [# particles of Z m < £ that hit B(0, R) c ] = P (£ m ' £ hits B(0, R) c ) , (3) 

where £ ro ' e denotes d-dimensional Brownian motion killed at rate e within T m and starting from 
under P(j. Known results about the survival of Brownian motion among soft obstacles (see in particular 
Sznitman's book |Szn98j) can be used to estimate the right-hand side of ©. On the other hand, one 
may expect that the inequality in ([3]) is "not too crude", but we do not pursue this matter here. 
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A major ingredient of the proof of Theorem [T]is the following quenched homogenization result. We 
need to introduce a rescaled version of the process Z ro,£ . For every e > and every t > 0, let us define 
a random measure X^' £ on ~R. d by setting, for every Borel subset A of M , 



xr ,e, 



z {A)=e J 2£%(dx)l e - U , A (x). (4) 

For every real x > 0, [x] denotes the integer part of x. 

Theorem 2. Except for a P '-negligible set of values ofw, the law of (X^ J ' £ )t>o under Pr e -i]$ converges 
weakly as e — > 0, in £/ie Skorokhod sense, to that of a super- Brownian motion with branching mechanism 
iP(k)( u ) := ir^ 2 + Kn started at 5o. 

The definition of super-Brownian motion with branching mechanism ipr K \ is recalled in Section 2 
below. 

As a hint of why Theorem [2] should be true, notice that for a given realization of the obstacles, the 
probability that a single Brownian motion starting from and killed at rate e within T m is still alive 
by time t > is given by 

E[exp{- e y l r „ (5) 

where £ denotes standard <i-dimensional Brownian motion. Let us focus on the integral within the 
exponential in ([5]). Averaging over the law of the obstacles and using Fubini's theorem, we obtain for 
each t > 



E 



E 



e [ lr w (&)da|l =eE[ / P[6 G I^ds 



EKt. 



We can thus guess, and easily prove, that the rescaled Brownian motion (i/e £ e -i^, i > 0), which 

is killed at rate 1 within y/e converges to Brownian motion killed at homogeneous rate k as 

s — > 0. Theorem [2] shows that an analogous convergence indeed holds in our more general framework 

of branching Brownian motions, for any fixed w contained in a set of P-probability one. 

Let us briefly explain how Theorem [1] is derived from Theorem |2j Consider a sequence (e n , R n ) such 

that R n — > oo and e n R\ converges to a positive constant o. By a simple scaling transformation, the 

l li 

probability P So (Z m > 6n hits R n A c ) coincides with F So (X™> 6n hits b n A c ), where b n = £n R n 

converges 

to \fa. We can then use Theorem [2] to investigate the asymptotic behaviour of the latter hitting prob- 
abilities. This limiting behaviour involves the corresponding hitting probabilities for super-Brownian 
motion, which are known to be related to solutions of semilinear partial differential equations from 
the work of Dynkin |Dy91, Dy93]. One difficulty in implementing the preceding idea comes from the 
fact that the convergence in Theorem |2] is not strong enough to ensure that hitting probabilities for 
the processes {X™ ,£ )t>o converge to hitting probabilities for the limiting process. Much of the proof 
of Theorem Q] in Section 4 is devoted to a precise justification of this property (Lemma |7J. 

To complete this introduction, let us mention that branching Brownian motion and superprocesses 
among random obstacles have been studied recently in several papers, including Englander and den 
Hollander [EdH03] and Englander |Eng08|. These papers concentrate on the case of supercritical 
branching, in contrast with critical branching which is considered here. See also the survey |Eng07| . A 
homogenization theorem related to Theorem [2] has been proved in [Veb09j for super-Brownian motion 
among hard obstacles, in the case when the intensity of the obstacles grows to infinity but their 
diameters shrink to 0. There is a huge literature about Brownian motion and random walks among 
(hard or soft) obstacles, and the reader may look at the book of Sznitman |Szn98] for additional 
references. 

The rest of this paper is laid out as follows. In Section (2j we introduce the basic notation and 
objects, and state several results about hitting probabilities for spatial branching processes we shall 
need in the sequel. Theorem [2] and Proposition [1] are proved in Section |3j Theorem [1] and Corollary [T] 
are then derived in Section HI 
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2 Preliminaries 



2.1 Notation 

We denote the set of all finite measures on M. d by M. t (R This set is equipped with the weak topology. 
We write _M p (K d ) for the subset of A4j(M. d ) which consists of all finite point measures on M. d . The 
integral of a function ip against the measure \x is denoted by (//, v 3 )- 

If E is a metric space, we denote the set of all bounded continuous functions on the space E by 
C(E) and we let ||/|| stand for the supremum norm of / G C(E). We write C 2 (M. d ) for the set of 
all twice continuously differentiable functions on R d , and C 2 (M. d ) for that of all bounded functions in 
C 2 (M. d ) whose first and second derivatives are also bounded. An index + added to this notation means 
that we require the functions to be nonnegative. We equip C 2 (M. d ) with the topology induced by the 
seminorms where for every R > 

\\f\\ {R) ■= sup !\f( x )\ + J2\^(x)\ + J2 l/rTf-to )• 

\x\<r{ ~t l9x i ' ^1 * J ' > 

If E is a Polish space, we let De[0, oo) be the set of all cadlag paths with values in E, equipped 
with the Skorokhod topology. We also let Ce[0, oo) stand for the set of all ^-valued continuous paths, 
which is equipped with the topology of uniform convergence on every compact set. 

If x G M. d and r > 0, B(x,r) denotes the open ball of radius r centered at x, and B(x,r) stands 
for the corresponding closed ball. More generally, the closure of a subset F of R^ is denoted by F. 
Lebesgue measure on R rf is denoted by A^- 

Finally, the notation £ = (S,t)t>o will stand for a standard Brownian motion in R rf , which starts 
from x under the probability measure P^. It will also be convenient to use the notation f^* 6 for 
Brownian motion in R d killed at rate e in the set As usual, the value of £ ro ' e after its killing time 
is a cemetery point A added to R d , and we agree that all functions vanish at A. 



2.2 Super- Brownian motion 

Let a > and set tp^(u) = (a 2 /2)u 2 + au, for every u > (the offspring distribution u, and thus the 
parameter a > are fixed throughout this work). Super-Brownian motion with branching mechanism 
■0(a) is the continuous strong Markov process with values in A^/(R rf ), whose transition kernels (Qt)t>0 
are characterized as follows: For every g G (7+(R a! ) and every [i G A4f(M. d ), we have for every t > 

Qt{n,dn') exp(- (//,#)) =exp(-(n,V t g)), (6) 

where the function Ut{x) = Vtg(x), t > 0, x G is the unique nonnegative solution of the semilinear 
parabolic problem 

— = \ Au - W) («) in (0, oo) x R d , 
^0 = 5- 

Let y = (Yt)f>o be a super-Brownian motion with branching mechanism V(a)' started at /U G 
^/(R '). Then, for every g G C^(R d ) 

e -aW _ e -0W _ J*(Y s ,-±Ag + ^ ia) (g)) e~^ds (7) 

is a martingale. It is well known that this martingale problem and the initial value /i characterize the 
law of Y. This is indeed an application of the classical "duality method" (see in particular Chapter 
4 in |EK 86j). The nonlinear semigroup g — > Vtg provides a deterministic dual to super-Brownian 
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motion, and the duality argument then shows that if a measure- valued process started from [i satisfies 
the preceding martingale problem, the Laplace functional of its value at time t must be given by the 
right-hand side of (j6|). See Section 1.6 of [EthOOj for more details. 

2.3 Branching Brownian motion among random obstacles 

In view of our applications (and in particular because we want to refer to some results of [Ch91|). 
it will be convenient to give a more formal description of the branching Brownian motions that were 
already introduced in Section [T] above. Recall that our offspring distribution v is assumed to be critical 
and that var(f) = a 2 G (0,oo). The probability generating function of v will be denoted by T. 

Let T be a Galton- Watson tree with offspring distribution v (see e.g. fLG05]). As usual, we view 
T as a random finite subset of 

oo 

U := (J N n , 

n=0 

where N = {1,2,3,...}, and N° = {0}. If v = (v 1 ,...,v n ) G U\{0}, the parent of v is denoted by 
v = («!,... ,v n -i) and we also use the notation v -< v' to mean that v' is a descendant of v distinct 
from v. Consider a collection (e v ,v G U) of independent exponential random variables with parameter 
1, which is also independent of T. We define for every v G U its birth time a v and its death time (3 V 
recursively by setting a = and (3 = e , and for every v G U\{0}, 

Q. v f3y , f3y (X V ~\- Gy . 

Let us now construct the spatial motions. Fix a starting point x G W 1 , and consider a collection 
(B v ,v G U) of independent standard Brownian motions in M. d (started from 0), independent of T and 
of (e v ,v G U). For every v G U, define the historical path uj v = (w",0 < t < f3 v ) associated with v 
in the following way. First uf = x + Bf for < t < /3 . Then, if v G U\{0}, set = ui^ for all 
< t < a v and 

<4 = <4 v +Bl-a v for a v <t< /3 V . 

A branching Brownian motion (without killing in obstacles) starting from 6 X is obtained by setting for 
every t > 0, 

Z t = 

where the notation v ~ t means that a v <t</3 v . 

In this formalism, it is now easy to introduce killing in obstacles. Consider yet another independent 
collection { r y v )veU of independent exponential random variables with parameter 1, and define for every 
w G f2, e > 0, and for every v G U 

C' £ := inf | s G [q„ p v ) : £ drl r „K) > e'S,}, 
where inf = oo. By setting 

%t ' £ = l{t<C^' £ and Cj' £ =oo, for every v'~<v} ^< , 

we obtain a branching Brownian motion killed at rate e in the obstacle set starting from 
An obvious extension of the preceding construction allows us to obtain branching Brownian motions 
starting from any point measure fi G Ai p (M. d ). 

We now recall a special case of the classical convergence of rescaled branching Brownian motions 
towards super-Brownian motion. For our applications, we state the case where particles are killed at a 
constant rate homogeneously over M. d (this case is obtained from the preceding construction of Z m,£ by 
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replacing T m by In the next two statements, for every e > 0, = [z\ )t>o denotes a branching 
Brownian motion with offspring distribution v, where particles are killed homogeneously over M. d at 
rate e. As previously Z^ starts from /i under the probability measure P^, for every /x S Ai p (R d ). 

Proposition 2. Let a > 0. For every e > 0, define a measure-valued process (X^ )^>o fry setting 



{XF,<p)=e J Z^> t (dx)^'x), 

for every cp G C(M d ). For every fixed rj > 0, f/ie lam of under P^-ij^ converges as e — > 0, in the 
Skorokhod sense, towards the law of super- Brownian motion with branching mechanism ip^ starting 
from r]5o. 

A proof of Proposition [2] can be found in Chapter 1 of [EthOO] in the case a = 0, and arguments 
are easily adapted to cover the general case. 

Finally, we shall use an estimate for the probability that a branching Brownian motion starting from 
5q exits a large ball centered at the origin. Similar estimates can be found in Sawyer and Fleischman 
|SF79| . but we provide a short proof for the sake of completeness. 

Lemma 1. Suppose that d = 1. There exist two positive constants C = Cq(v) and C[ = C{(^) such 
that, for every e 6 [0, 1] and r > 1, 

C',(r- 2 l {r ^ } +ee- r ^l {r> ^ } ) < P*(Z« hits (-r,r) c ) < C[ (V 2 l {r ^ } + ee~ r ^ l {r> ^ } 

Remark. In dimension d, we can apply the bounds of Lemma [T] to the projection of Z^ on each 
coordinate. In particular, taking e = 0, we deduce from the upper bound of the lemma that we have 
for every r > 0, 

F So hits B(0, r) c ) < C'{ (r + 1)~ 2 (8) 

with a constant C'{ = C'{{d,v). 

Proof. It clearly suffices to prove that the stated bounds hold for the quantity ¥g r (Z^ hits (— oo,0]) 
instead of ¥s (Z^ hits (— r, r) c ). We fix e > 0, and for every x > and t > 0, we set 

q e (x,t) = F$ x {z^ does not hit (— oo,0] before time t 

and 

p e {x) = F Sx (z& hits (-oo,0]) =limt (l-q E (x,t)). 

In this proof only, we write P| for the probability under which £ is a Brownian motion starting from 
x and killed at rate e (upon killing, £ is sent to the cemetery point A and we recall that all functions 
vanish at A). Write S := inf{i > : £t £ (— oo,0]}. By standard arguments (see e.g. the proof of 
Proposition II. 3 in |LG99| ). the function q e solves the integral equation 

q e (x,t)=P x (S>t)+E x / (r(q E (Z s ,t-s))-q £ (ti s ,t-s))ds 

JO 

where we recall that T denotes the generating function of the offspring distribution v. For every 
a S [0, 1], set $(a) = T(l — a) — (1 — a). Note that $(0) = and the function 3> is monotone increasing 
under our assumptions. Furthermore, <5(a) = (o" 2 /2)a 2 + o(a 2 ) when a — > 0. By a monotone passage 
to the limit we get that, for every x > 0, 

p £ (x) + E x \ f 3(p e (&)) ds] = P X (S < oo). (9) 
Jo J 
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From standard connections between Brownian motion and partial differential equations (see e.g. Chap- 
ter 8 of Durrett [Dur84j), it follows that the function p £ satisfies the differential equation 

\p'e =ep £ + $(p £ ) 

on (0,oo) with boundary conditions p e (0) = l,p e (oo) = 0. By solving this differential equation, we 
get, for every x > 0, 

du 

- x 



lp e ( x ) yj2eu 2 + 4T(u) 

where T(u) = J " ) dv . Note that there exist positive constants c, d such that cu 3 < T(u) < c'u 5 
for every u G [0, 1]. The desired bounds then follow from easy analytic arguments. □ 



2.4 Hitting probabilities for super-Brownian motion 

Let Y"(°) = (Y t {a > )t>o be a super-Brownian motion with branching mechanism ip^ for some a > 0. 
Suppose that starts from 11 under the probability measure P^, for every 11 S Mf(R d ). 
Following [Per 02 1 p. 200], we define the range of y» by 



K(Y^) = \J Mj su PP(^ (a) ) > 

e>0 \t=e J 

where for every /i £ A4j(M. d ), supp(/x) denotes the topological support of fx. 

Let D be a domain in M. d and let x £ D. Consider the process Y^ started from 5 X . We say 
that y( a ) hits D c if the range TZ(Y^) intersects D c . By a famous result of Dynkin |Dy91| [Dy93| the 
function 

uf a) (x) = -log (l - PsAny (a) ) n D c ± 0)) , x € D, 

is the maximal nonnegative solution of the semilinear partial differential equation ^Aii = if)^(u) in 
D. 

Under mild regularity assumptions on D (which hold e.g. when D satisfies an exterior cone condition 
at every point of dD), the function vP-, has boundary value +00 at every point of dD and is the unique 

nonnegative solution of the equation ^Ati = if)^(u) in D with boundary value +00 everywhere on 
dD. A discussion of this result and related ones can be found in Chapter VI of the book [LG99J. This 
reference considers only the case a = 0, but the same results can be obtained for any a > by similar 
arguments: Note that the Brownian snake approach can be extended from the case a = considered in 
[LG99J to a > 0, simply by replacing the reflecting Brownian motion driving the snake by a reflecting 
Brownian motion with negative drift (see Chapter 4 of |DLG02J for a discussion of the snake approach 
to superprocesses with a general branching mechanism). 

We shall be interested in the special case D = A. Recall that 6 A and that we assume A is a 
domain of class C 2 , meaning that the boundary of A can be represented locally as the graph of a twice 
continuously differentiable function, in a suitable system of coordinates. We write u^(x) = u^ a Ax) to 
simplify notation. From the analytic viewpoint, the function ut a \ may be constructed as follows. For 
every integer n > 1, let Ui a \ n be the unique nonnegative solution of the nonlinear Dirichlet problem 

f \Au = i/}(a)(u) in A, 
1 u\dA = n . 

Then uu) = lim f Uf a \ n as n — > 00. 

The following lemma records certain analytic properties which will be useful in the forthcoming 
proofs. 
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Lemma 2. (i) Let x G A. The function a — > ur a \(x) is continuous and nonincreasing on [0,oo), and 
tends to as a —> oo. 

(ii) For every 5 G (0, dist(0, A c )), let A$ be the subdomain of A defined as the connected component of 
the open set {x G A : dist(x,A c ) > 5} that contains 0. Then, for every a > 0, it/ \(0) tends to U( a )(0) 
as 5 — > 0. 

Proof, (i) Let us first verify that the function a — > uu)(x) is monotone nonincreasing, for every 
x G A. To see this, we apply a standard comparison principle (see e.g. Lemma V.7 in [LG99]) to 
obtain that Ur a i\ n < Ur a \ n if a < a' , for every n > 1. It then suffices to let n — >• oo. 

Let (a/c)fc>i be a sequence of nonnegative reals increasing to a G (0, oo). We can set for every x £ A 

v(x) = lim | ur ak \(x), 

and we have v > U( ). In order to verify that v < ii( ), we only need to check that v solves ^Av = -0(a) ( v ) 
in A (recall that ut a \ is the maximal nonnegative solution of this equation). To do so, let B be an 
open ball whose closure B is contained in A. For every k > 1, the restriction of ut a ,) to B solves the 
equation ^Au = V'(a fc )( n ) i n B- By the probabilistic interpretation of the integral equation associated 
with this PDE (see e.g. Chapter V in [LG99]), this implies that, for every x G B, 



u 



(flk 



)(x) + E a 







lp(a k )(u { a k ){£s))ds 



where we recall our notation ^ for a Brownian motion starting from x under the probability measure 
P x , and tb '■= inf{i > : & ^ -B}. By passing to the limit fc — >• oo in the previous display, we can 
write 



v(x) + E x 







which is enough to obtain that v solves \Av = V ; (a)(' u ) m B, and therefore in A since B was arbitrary. 

Similar arguments show that, if (afc)fc>i is a decreasing sequence of nonnegative reals converging 
to a G [0, oo), then ur ak )(x) converges to uu\(x) for every x G A. Finally, the fact that ur a \(x) tends 
to as a —t oo can be obtained from the comparison principle: If B is a ball such that B C A, the 
restriction of ui a \ to B is bounded above by the solution vi a \ of the linear equation ^Av( a ) = ao( a ) 
in with boundary value equal to the restriction of um\ to B. It is easily seen that v/ a \(x) — > as 
a — > oo, for instance by using the Feynman-Kac formula. 

(ii) Fix a > 0. If < 5 < 5', the closure of A^/ is contained in Ag. The restriction of u,J[ to ^4^/ is 

a nonnegative solution of ^Ait = ip^(u) in Aj/ and is thus bounded above by . Hence, for every 

fixed x G A the function 5 —r u,{(x), which is defined for 5 > small enough, is nondecreasing and 
we can set 

v(x) = \imu^\(x). 

By the same argument we used to obtain the monotonicity of the mapping 5 — > u^A(x), we also have 
v(x) > U( a )(x) for every x G A. To obtain the reverse inequality v < Ur a \, it is enough to verify that 
v solves \Av = t/jr a )(v) in B. But this follows by arguments similar to those we used in the proof of 
part (i) of the lemma. □ 

Lemma 3. For every a > and x G A, 

{K(Y^) n A c ^ 0} = {K(Y^) n {Af ^ 0} , P 5x a.s. 
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Proof. The inclusion 

{7£(Y (a) ) nA c ^(D}D {^(Y (a) ) n (A) c + 0} 

is trivial. To show the reverse inclusion, we may argue as follows. By Theorem IV. 9 in [LG99J 
(which holds under much less stringent assumptions on A), the event {TZ(Y^) D A c ^ 0} holds if and 
only if the exit measure of the super-Brownian motion from A is nonzero. Applying the special 
Markov property of superprocesses |Dy93 Theorem 1.3], we see that it is enough to prove that for 



super-Brownian motion starting from a nonzero initial measure supported on dA, the support process 
immediately hits (A) c . Under our regularity assumptions on A, this easily follows from Wiener's test 
for super-Brownian motion (Theorem VI. 6 in [LG99]). □ 

From now on, we write {y( a ) hits F} for the event {TZ(Y^) n F / 0}. 

3 Quenched convergence to super-Brownian motion 

The main goal of this section is to prove Theorem [2] At the end of the section, we also establish 
Proposition [TJ using certain arguments related to the proof of Theorem [2] To simplify notation, we 
set for every w G ft and e G (0, 1), 

iei 

The following lemma identifies a martingale problem solved by our branching Brownian motion 
Z m ' £ . It can be proved by standard arguments: Since killing in the obstacles can be viewed as a 
location-dependent branching phenomenon, our martingale problem is indeed the same as in (4.10) of 
[EK861 Section 9.4] (see also Section 1.2 in [EthOO] ). 

Lemma 4. Let w G f2 and e > 0. Under each probability P„ ; fi G A4 p (M. d ), the process Z ro,£ solves 
the following martingale problem: For every f G C+(M. d ) such that < inf / < / < 1, the process 



is a martingale. 

We can derive from Lemma U (or from a direct argument) that for every g G C 2 (M. d ), 

M t (g) := (Z?*,g) ~ (Zo' £ ,g) - f {zf> £ , Ug - el r ,g) ds (10) 
is a martingale. An easy computation gives that the square bracket of this martingale is 

[M(g),M(g)] t = [\z^,Vg.Vg) ds + £ (Z?* - Z?_! £ ,g) 2 . (11) 

J ° 0<s<t 

The last sum in the right-hand side is an increasing process with compensator 

[\z?> £ ,(o- 2 + el r Jg 2 )ds. (12) 
J o 

The proof of Theorem [2] relies on the following two results, in which we use the notation (3^t)t>o f° r 
the canonical process on Dj^ [0, oo). Recall from (|4]) the definition of the process X m ' £ in terms 
ofZ ro < £ . 
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Lemma 5. For every w G ft and e G (0,1), let n ro,e 6e i/ie law of the process X ro,e under ii<j • 
(i) For ewery 5 > and T > 0, there is a compact subset ofM. d such that, for every w G CI, 



sup n ro ' e sup y t (^ T ) < s 

££(0,1) v 0<t<T ' y 

(ii) For every g G C ,2 (M d ) andw G f2, £/ie collection of the laws of the process ((3^,<?))i>o under n ro ' e , 
e G (0, 1), is relatively compact in the space of all probability measures on D%[0,oo). 

Consequently, for every w G O, i/ie collection (n ro,£ ) eg ( ,i) * s "relatively compact in the space of all 
probability measures on Dj^ .(Rd) [0, oo) . 

The last assertion of the lemma is an immediate consequence of (i) and (ii) using Theorem II. 4.1 
in |Per02j . 

Proposition 3. Let g G C?(R ). There exists a measurable subset Sl g of Q such that P(f2 g ) = 1 and 
i/ie following holds for every w G f2 s . For ewery s,i > 0, for every integer p G N and ewery choice of 
tx,...,t p G [0,i] and fi, - ■ ■ ,f p G (^(.M/QR^)), we acwe 



limE [e -i ]( $ 



We postpone the proof of Lemma |5] and Proposition [3j and explain how Theorem [2] follows from 
these two statements. We choose a countable dense subset G of C+(M. d ) and set 

n' = n *V 

Fix ti7 G fi'. By Lemma [5j the collection (n ro,e ) £g (o i i) is relatively compact. Let IT* be a sequential 
limit of this collection as e tends to 0. We deduce from Proposition [3] that, for every g G G, for every 
s, t > and every choice of t 1} . . . , t p G [0, t] and / x , . . . , f p G C(X /(R d )) , we have 



IT 



(34, -\&g + i>( K) (g)) e-W^du} J[ 



(13) 



The required passage to the limit under the expectation sign is easily justified: Note that, by comparing 
with the case when there is no killing and using standard results of the theory of critical branching 
processes, we have for every T > 

sup (E [e -i 5o] [ sup (X?' £ ,l) 2 ]) <oo. (14) 

ee(0,l) V L 0<r<T J/ 

Since G is dense in C^_(M. d ), another easily justified passage to the limit shows that f 1 1 3 p holds for every 
g G C+(M. d ). Thus, IT* satisfies the martingale problem for super-Brownian motion with branching 
mechanism as stated in Section |2] Since it is also clear that IT*(Yo = 5q) = 1, IT* must be the law 
of super-Brownian motion with branching mechanism ipr K \ started from 5o. This completes the proof 
of Theorem [2j but we still need to establish Lemma [5] and Proposition |3] □ 

Proof of Lemma [5l The compact containment condition (i) in the lemma is immediately obtained 
by observing that X™' £ is dominated by X° = X m '° and by using the convergence of rescaled branching 
Brownian motions (without killing in the obstacles) towards super-Brownian motion. So we only need 
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to verify (ii). In the remaining part of the proof, we fix w G fi and g G Cl(R ). To simplify notation, 
we set 

Xl = {X^,g). 

By the remarks following Lemma [4] and an elementary scaling transformation, we have 

XZ = X% + Mt + V t E , 

where 

Vf = j* (Xft^g-lv^g) ds 
and M e is a martingale, whose square bracket is given by 

[M £ ,M% = f(XT' £ ,eVg.Vg) ds + £ (X?> s - X^ £ ,g) 2 . 

J° 0<s<t 

Furthermore, the oblique bracket of M e is equal to 

(M £ ,M £ ) t = I (Xf-VV + e{Vg.Vg + l r ^ 2 )} ds. 
Jo 

By standard criteria (see in particular Theorem VI. 4. 13 in Jacod and Shiryaev [ JS87J ) , the tightness 
of the laws of the processes X s , e G (0, 1), will follow if we can verify that the laws of the processes V s 
and (M e ,M e ), for e G (0, 1), are tight in the space of probability measures on Cr[0, oo). But this is 
immediate from the preceding explicit formulas and fjl4[) . □ 

Proof of Proposition [H Let us fix s, t > 0, h, . . . , t p G [0,t] and fi,...,f p G C(M f (M d )). Let 
vj G Q and e G (0, 1). By Lemma |4] applied to the function f(x) = exp j — eg{xy/e)\ , we can write 



E 



'x^,^(-e 2 Ag + e 3 Vg.Vg) 



+ el r ^{l - e- £9 )e £9 + e e9 (T(e- £9 ) - e^ 9 )) e 



V 



E 



* +s ' I.e. 



X^,--Ag + -Vg.Vg 



' i=l 



(15) 



where r/ E := £ _1 (1 - e" £9 )e £9 - 5 and fj e := (2/a 2 )e~ 2 e S9 (T(e- £9 ) - e' 69 ) - g 2 . Using a Taylor 
expansion and the facts that T(l) = T'(l) = 1, T"(l) = a 2 , we readily obtain that for e small 
enough we have \\rj £ \\ < C\{g)e for some constant C\{g) depending only on g, and ||ry e || — > as e — > 0. 
Now, using Fubini's theorem, the fact that ((X^' £ , l)) r >o is a super martingale and the inequality 
E [e -i ]4b [(X^ £ , 1)] = eie- 1 } < 1, we have 



E 



J i=l 

< (|(liv 5 .v 5 || + 2d( 5 )) + ff||*y) (f[ H/iii^V 1 !* [/ +S(X "' £ ' 1 

< (C 2 (g) e + (a 2 /2)\\U) (f[\\M\)s 



du 



i=l 
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with a constant C%(g) depending only on g. On the other hand, 



E 



[ V- la^e-^dujn/*^) 
• 7t J 1=1 

J t t+S (X^,K 9 ) e-^'^du} n/^XH 



(16) 



The absolute value of the second term in the right-hand side of (fl6|) is bounded above by 



n 

i=i 



t+s 



E [e- 1 ]<5o 



fX ro ' £ ,flT 



du. 



Consequently, going back to fj 1 5 [) and using the preceding estimates, we obtain for every zu £ CI and 
£€(0,1), 



E [e- 1 ]5o 



8=1 



< (nil/ill) ({C , 2(g)e + (ff 2 /2)||%[|}5 + r e (tu,(7,t ) t + S ) 



(17) 



where r e (zu, g,t\,t2) is defined for < ii < ty, by 

-t 2 



Lemma 6. Let u > 0. TTien, P-a.s. 



ti 



1™ - K)g) 



du. 



X>*, (lr.-K)g) 



(18) 



Assuming that the lemma is proved, we can readily complete the proof of Proposition [3] Using 
Fubini's theorem, we may find a set Cl g with P[f2 ff ] = 1, such that, for every w £ Cl g , the convergence 
in (|18p holds simultaneously for all u > 0, except possibly on a set of zero Lebesgue measure (depending 
on w). Since for every w £ CI, e £ (0, 1) and u > 



E 



\(X^, {1^-k) 9 )\ <|H|E [£ - 1](5o [(X^ £ ,l)]<|| g || 



dominated convergence guarantees that for each w £ Cl g , lim^o r e( ro ) <7) t } t + s) =0 for all t, s > 0. It 
follows that the right-hand side of (|17p tends to as e — > when w 6 fi s , which completes the proof 
of Proposition [3] □ 

Proof of Lemma O Recall that £ ro ' e denotes standard (i-dimensional Brownian motion killed at 
rate e within r ro . We also use the notation 101 Brownian motion killed at rate 1 within T^. Both 
processes and x ro ' e start from x under the probability measure P^. 

We first recall classical moment formulas for branching Brownian motion. For every x £ M. d , k £ N, 
and every bounded measurable function h on R d , we have 



E fc u<zr^>] = ^[Mr* £ )] 
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and 



E k s x [(Z?' £ ,hy]=kE x [h(Cn 2 ] + Hk-l)[VM&')]) +ka 2 E 



'(fa (V £T ,e[h(Z?4)] S2 



These formulas are easily derived, first for k = 1, from the well-known formula for the Laplace functional 
of Z™' £ under P^: See e.g. Proposition II. 3 in [LG99J. 

Recalling the definition of X m ' £ in terms of Z™' 6 , we get similar formulas for the first and second 
moment of {Xf' £ , h). In particular, for s > and for every w G fi and e £ (0,1), we have 



E 



W' £ , (li» - K )g) 



( [e -i]2 _ ^ (0 , S) 5) 2 + e[e -i]^ ( o, s , 5 ), (19) 



where, for every x G 



^ £ (x, S , 5 ) := E4(lr^(xr ,£ ) " «0<7(x7' e )] : 



and 



A"' £ (x,s,5() := cr 2 E a 



+ eE4(lr^(xr e )-^)Vx 



We claim that, for every x G M. d and s > 0, 



P— a.s., ]imAT' e (x,s,g)=0. 



(20) 



Assume for the moment that the claim holds. By taking x = in (|20p . we obtain that A™' £ (0, s,g) 
tends to as e — > 0, P-a.s. Consider next A^' £ (0, s, g). The second term in the formula for A^' £ (x, s, g) 
obviously tends to uniformly in x, s and independently of the obstacles, as e — > 0. To handle the 
first term, use Fubini's theorem to obtain that, for every w belonging to a set ft of full probability, 
there exists a set A/" ro C M. d x M + of zero Lebesgue measure such that the convergence in (I20p holds 
simultaneously for all (x, s) ^ M m . We can then write the first term in the formula for A^' e (0, s,g) as 
follows: 



rr~ I dv 

'o 



dy vT £ {^v) V y 



, vXs 



K)g(x7-v) 



(21) 



where p^' £ (-, •) is the transition kernel of x ro,£ a t time v. Plainly, we have p^' £ (x, y) < p v (x, y) for each 
v > and e G (0,1), where p v (-,-) is the transition kernel of standard Brownian motion at time v. 
Using this bound and dominated convergence gives us that the quantity in (|2ip tends to as e — > 0, 
for every w G ft. Hence A™' e (Q,s,g) tends to as e — > 0, P-a.s. The result of Lemma |6] now follows 
from (|19p and the Cauchy-Schwarz inequality. 

We still have to prove our claim (l20|) . We fix x G M rf and s > 0. Let w G fi and e G (0, 1). From 
the definition of Brownian motion killed at rate 1 in rfL, we have 



E, 



K)g(x 



[lr«, (6) ~ ex P { - J ^(Z^duj 



(22) 



Let r\ > and choose G (0, s) such that e 9 > 1 — r\. By the Markov property applied at time s — 6, 
the right-hand side of (I22p is equal to 



E, 



exp 

+E,; 



Ire 



E 



exp 



exp 



lr| 

-0 

lr 



(£ u )du\ E 



lr* (6) - k)9(£o) exp { - ^ l rfo (£„)cfo} 

(li* (6) -«)»(&)] 
I- . IdrcJZe) ~ *)sfe)(exp { - f l a (£«)£*«} 



(23) 
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The condition 1 — e 8 < r\ entails that the absolute value of the second term in the right-hand side of 
(|23p is bounded above by \\g\\f]. Suppose we know that for every y G M. d we have 



P— a.s., 



lim E 

e^O 



0. 



(24) 



Then, using the fact that the law of £ s _# is absolutely continuous with respect to Lebesgue measure 
(and Fubini's theorem to see that the convergence in (|24p holds simultaneously for almost all y G M. d , 
P-a.s.), we conclude that the first term in the right-hand side of (123[) converges to with P-probability 
1 as e — > 0. Consequently, we have P-a.s. 



limsup E* {Irejx7' £ )-K)9(x7' e ) 

£— >-0 L 



< \\g\\v- 



Since rj was arbitrary, our claim (|20p follows. 

It thus remains to prove that (124p holds. We fix y G M. d and 9 > 0. In the following, £' stands for 
another Brownian motion independent of £, which also starts from y under the probability measure 
~P y . For each e G (0, 1), we have 



E 



E v [(lr«,(&) 



E 

E, 



E,[(lr^(6) - (&) - k)<7(&)</(# 



(25) 



where the last line uses Fubini's theorem and the definition of k. Recall that the measure O is 
supported on compact sets which are contained in the fixed ball B(0, ro). If \x — x'\ > 2ro, the sets 
{(z, K) : x G z + K and K C B(0, r )} and {(z, K) : x' G z + K and K C B(0, r )} are disjoint, and 
so the events {x G and {V G are independent under P. Recalling that = -y/e we see 
that if \x — x'\ > 2rQyfe we have P[x G T^; x' G T^] = P[x G P^jPfa;' G T^] = k 2 , which enables us 
to write 



Going back to 025 p . we obtain 



E 



E,[(i^(e e )-«)<?(^)] : 



: E,, 



<N| 2 Py[l&-&l <2r v^] 
= ||5l| 2 Po[|e 2 e| < 2r 0V ^] 



.d/2 



(26) 



where the constant C > depends only on rg. Let > 0. By the Markov inequality, we can write 



|E„[(lr-,(6>) - I > n 



< C 



T d/2 



Applying the last bound with e = = k 3 , for every k G N, yields a convergent series, and so by the 
Borel-Cantelli lemma we have P-a.s., 



lim sup 

k— >oo 



Ey[(W&)-«)<7(&) 



< 



Since rj > was arbitrary, the convergence (|24p holds along the sequence (£fc)fc>i- 
To complete the proof, we set for every w G fl and e G (0, 1) 
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(recall that y and 9 are fixed). We shall prove the following result: For every zu £ ft, 

lim sup | U m ^ e - U^ )Bk | = 0. (27) 

k ^°° (fc+l)-3<e<fc-3 

Combining (|27p with the fact that the convergence (|24p is true along the sequence (£fc)fc>i will then 
lead to the desired conclusion. In order to prove (|27p . we first note that for every e G (0, 1), 

1 f \*-y\ 2 f 

U VJE = - — tjtt / Ire (x) q(x)e 292 d x — I l T e (x)h(x)dx, 

with a function h E C 2 _{M. d ) which depends on y and 9, but not on w. Furthermore, for any fixed 
r] > we can find a large closed ball B centered at the origin and such that J BC h(x)dx < rj. Hence, if 
we set 

U^ £ = [ h(x) l r ,Jx)dx, 

we have \U Wie — U'^ £ \ < r\ for every e E (0, 1). Thanks to this remark, it is enough to prove that (12?]) 
holds when U m%e and U^ }£k are replaced by e and E respectively. 
Let k E N and s G [(k + 1)~ 3 , A; -3 ). We have 

f4, E =e d/2 / h{x^)l T ^x)dx 
Je- 1 / 2 B 

= e d/2 [ h(x^)l T ^{x)dx + e d/2 [ h(x^)l r ^(x)dx. (28) 

Jk 3 / 2 B J(e- 1 / 2 B)\(k 3 / 2 B) 
Now, the first term in the right-hand side of (|28p is equal to 

where 

i(zu,e,k) = (ek 3 ) d / 2 [ \h((k 3 / 2 V~e)y) - h(y)] l v ^' 2 y)dy. 

JB 1 i 

Note that < 1 — k 3 l 2 \fe < £ with a constant C independent of e and fc, from which it easily follows 
that 

c 

sup \c(w,e,k)\ < — \\Vh\\, 
(fc+i)- 3 <e<fc- 3 k 

with a constant C" depending only on B. 

Similarly, the second term in the right-hand side of (|28p is bounded above by 

e d l 2 \\h\\[e~ d ^ 2 - k zd l 2 ]\ d {B) < C " 



k 

Finally, from (|28p and the preceding estimates, we have 

i^-^j< (1 - (£A¥ / V ^ + ^ + qM, 

and the convergence (i27|) follows (note that the sequence (U^ £ )k>i is bounded by construction). This 
completes the proof. □ 

Proof of Proposition [TJ Clearly we get a lower bound for the quantity ¥^ {Z w,e hits -B(0, R) c ) if we 
replace Z ro ' e by branching Brownian motion killed homogeneously over M. d at rate e. This intuitively 
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obvious fact can be proved by a coupling argument. Part (i) of the proposition thus follows by applying 
the lower bound in Lemma [1] to the projection of on the first coordinate. 

Let us turn to the proof of the first assertion in (ii). We may restrict our attention to the values 
of R > 1 and e £ [0, 1] such that E\f~R is larger than a fixed constant depending on d: Just observe 
that the probability that (unkilled) branching Brownian motion starting with one initial particle at 
the origin exits the unit ball is bounded above by a constant (depending on v) strictly smaller than 1. 

We then have the following inequality: For every w G Q, and e £ (0, 1), 



5o [Z m > e hits B(0, R) c ] < P [C £ hits B(0, R) c ] . 



(29) 



Indeed, using the formalism of Subsection 12. 3| the criticality of the offspring distribution can be used 
to check that the right-hand side of (|29p is just the expected value of the number of those historical 
paths uj v that first exit B(0,R) during the interval [a v ,/3 v A (1?' £ ). Alternatively, it is easy to derive 
an integral equation similar to (j9|) for the function x — > ¥s x [Z m,£ hits i?(0,i?) c ], and the bound (|29|) 
then trivially follows from this integral equation. 

Then, let us bound P ® P [£ ro ' e hits B(0, R) c ] . We write £ ro ' e ' fc for the fc-th coordinate of £ ro ' e , for 
every k = 1, . . . , d. We observe that 



F [C' £ hits B(0,R) c ] < ^ P <g> P [C' e,k hits (-R/Vd,R/Vd) c 



(30) 



fc=i 



Let us consider the first term in the sum. Define N £ r := [Ry/e/yd] and for every j £ Z, set 
T? = inf{i > : ^ = je^ 1 ^ 2 }, where £j stands for the first coordinate of With this notation, we 
can write 

P® Poff 7 ' 2 * 1 hits (-R/Vd,R/Vd) c ] 

r^' 1 hits(-iY £ ,^- 1 /2 )A r £i?e -i/2 )C 



<P® P 
<E®E 



exp 



+ E(g>En 



exp 



(31) 



Consider the first term in the right-hand side of (I3ip . Recall the definition ([I]) of the set of obstacles, 
and set for every j £ N and e > 0, 



td-l 



e) := (J (a* + Xi) and T^(j) := ^ r ro (j, e). 

Note that the random sets e), j £ N, are independent under P, by properties of Poisson measures. 

We have then 



E® En 



exp | 









<E(g)E 



:E(g)E n 



exp 



exp 



=1 



j=i ^^-i 



lr„y, e )(C«)rf« 



L r ro (i, £ ) 



(32) 



where the equality comes from an application of the strong Markov property of £, together with the 
independence of the random sets r ro (j, e) and the fact that the distribution of each of these random 
sets is invariant under translations by elements of {0} x By scaling, if T\ denotes the entrance 

time of £ into [l,oo) x we can write 

E(g>E expl-ef lr ro (i )£ )(6)ds) = E(g>E exp j - f l v ^ {1) {i s )ds 



a. 
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We then observe that 



Ti 



ii*(i) (6)i(o,i) (tl) ds — 1 



£->0 



l(o,i) 



(33) 



where the notation — > refers to convergence in probability under P <g> Po (and we recall that is the 

first coordinate of £ s ). To see this, we use arguments similar to the proof of Lemma |6] Notice that 
P[y G r^(l)] < K if y G (0, 1) X W^ 1 , with equality if y G (roy 7 ^ 1 — tq^/E) x R d_1 . Using this remark, 
and the same argument as in the derivation of (I26p . we can write for every fixed u > 0, 



E<g>E 



T-lAu 



: E (g) En 



rliAU rl\AU 

/ / (li* (i)(6) " «)(lre,(i)(6) - «) 1(0,1) Ki) 1(0,1) fe 1 ) dsdt 
Jo Jo 



— Eo 
= Eo 
<E 



TiAn pTiAu 









rTiAu rTiAu ^ ^ ^ _ ^ &1) 

JO ■/ 

+ 0(e 1 / 2 ) 



{P[6,6 G r^(l)] - k 2 } l m (e s )l(o,i)($)dsdt 



+ 0(e 1 ' 2 ) 



dsdt 



+ 0(e 



l/2\ 



Jo 



1 {|6-?t|<2rov / ^} dsdt 



where the error term 0{e 1 ^ 2 ) corresponds to the contribution of times s,t such that ^ or belongs 
to the set (0,roy/e\ U [1 — roy/E, 1). The preceding quantity tends to as e — > 0, which yields the 
convergence (|33p . Since the limiting variable in (j33[) is (strictly) positive a.s., we can find S\ G (0, 1) 
and c\ < 1 such that a e < c\ for every e G (0, ex)' Using (13 ip and (132 p . we see that the first term 
of the sum in (|30p is bounded above by 2c x e,R provided that e G (0,£i). The same bound will hold, 
possibly with a different value of e 1; for the other terms of the sum in (|30p . Therefore, we arrive at 
the bound 

P ® Po[£ ro ' e hits B(0,R) c ] < 2dcf £ ' fl , 

for e G (0,£i). Now recall from the beginning of the proof that we may restrict our attention to the 
values of R and e such that £y/R is larger than a fixed constant depending on d. The first assertion in 
(ii) now follows from (12 9 p and the last bound. 

Let us turn to the second assertion. Let C 2 > be a positive constant whose choice will be specified 
later. By simple comparison arguments, it is enough to prove the desired estimate when R is of the 
form R = 2 k , for k G N large enough, and e is of the form e = 2 _ - ? , with j G {0, 1, . . . , 2k} such that 

£>C 2 (loglogi?) 2 AR 2 . 

By the first assertion in (ii) and the Markov inequality, 

F So [Z™' £ hits 5(0, R) c ] > exp(-Cii? v / e/2)j < exp(-CiRy/e/2). (34) 

However, if R = 2 k and e > C 2 (loglog J R) 2 /i? 2 , we have 

R\/E > \fC~2 log log R = v 7 ^ (log k + log log 2) . 

Using this bound, we can choose the constant C 2 sufficiently large so that we get a convergent series 
when we sum the right-hand side of (|34[) over all R = 2 k and e = 2 _ - ? for j G {0, 1, . . . , 2k} such that 
£ > C 2 (loglog R) 2 j R 2 . The Borel-Cantelli lemma now yields the desired result. 

□ 
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4 Proof of the main result 



In this section, we prove Theorem [TJ We fix the environment w such that the weak convergence of 
Theorem [2] holds, and derive the convergence in Theorem [1] for this fixed value of the environment. For 
the sake of simplicity, we shall omit w in the notation and write Z £ instead of and X s instead 

of X ro ' £ . 

We shall verify that for any increasing sequence (R n ) n >i °f positive reals converging to +oo and 
any sequence (e n ) n >i of nonnegative reals such that e n R\ — >• a G [0, oo], we have 



lim R 2 n F So {Z £ " hits R n A c ) = u (Ka) (0), (35) 



where ^(oo)(0) = by convention. 

The statement of Theorem [T] follows from this convergence. Indeed, if the conclusion of the theorem 
fails, then we can find a sequence R n f oo and a sequence (e ra ) of nonnegative reals such that, for every 
n > 1, 

\R 2 n F 5o (Z £ " hits R n A c ) - u (KenRl) (0)\ > 5 

for some constant 5 > 0. By extracting a subsequence, we may assume that £ n Rn — > o G [0, oo] and 
thus obtain a contradiction with (|35p since we know from Lemma [2] that the mapping b — > Un,) (0) is 
continuous on [0,oo]. 

In proving (|35p . we may assume that e n — > as n — > oo. Indeed, suppose that fj35 1) holds in this 
particular case and let (e' n ) be a sequence that does not converge to 0. If the sequence e^R^ converges 
then necessarily its limit is +oo, and we can find another sequence such that < e'^ < e' n , e" — > 
and e" n R\ — > oo. So, if we know that (j35j) holds in the case when the sequence (e n ) tends to 0, we 
obtain 

lim RlF So (Z £ » hits R n A c ) = 0. 



However, from the inequality < e' n and a coupling argument (obvious if one uses the construction 
described in Subsection 12. 3ft . we get the same result for the sequence (s' n ). 

A similar comparison argument shows that it is enough to prove (|35[) in the case when a < oo. 
Otherwise, it suffices to replace e n by e n A bR~ 2 and let b — > oo, using the fact that ujJO) — > as 
b — > oo. 

Let us now proceed to the proof of 035 p . We fix the sequences R n f oo and e n — > such that 
e n R\ — > a G [0, oo). We first assume that a > 0. The case a = will be discussed at the end of the 
section. 

Let B be a closed subset of M. d . For every e > 0, we have by the definition of X s 

F [£ - 1]5o (X £ hits B) = ¥ [E - 1]8o (3t > : Xf(B) > 0) 

= F [£ -i ]So (3t>0: J Z e £ „ H {dx) > 

= P [e - 1]5o (Z £ hits e-^B) 
= 1 - F So {Z £ does not hit e _1 / 2 5)[ e_1 ] ; 



since (for a fixed environment) the law of Z e under P[ e -i],j is obtained by adding [e x ] independent 

1/2 

copies of Z e under F$ . Applying the preceding identity with e = e n and B = b n A c , where b n = e„ R n , 
gives us that 

1 - P So (Z £n does not hit R n A c )^ = P [£ -i ]So (X 6n hits b n A c ) (36) 

By Theorem [2J we know that the law of X 6n under P^-ij^ converges as n — > oo to the law of super- 
Brownian motion with branching mechanism ipr K ^ started from 5q. The next lemma is essentially a 
consequence of this convergence. We use the notation of Subsection 12.41 
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Lemma 7. We have 

lim P, £ -i ]s (X £ » hits b n A c ) = P 5o (F (K) hits bA c ), 

where b = y/a = lim 6 n . 

We postpone the proof of Lemma [7] and proceed to the proof of (|35p . in the case when a > 0. By 
the results recalled in Subsection 12.41 we know that 



P 5o (y (K) hits bA c ) = 1 - exp(-u(0)) , (37) 

where the function (v(x),x £ bA) is the unique nonnegative solution of the singular boundary value 
problem 

JiAv = V( K )(u) in bA , 

\ U\d{bA) = +oo . 

It is immediate to verify that U( Ka -j(x) = av(bx) for every x 6 A, and in particular ur Ka \(0) = av(0). 
From (|36p . (|37p and Lemma [JJ we obtain 

lim (1 - F So {Z £ " hits R n A c ))^ = exp(-v(0)) 

and thus 



lim E^FgJZ^ hits R n A c ) = v(0), 

n - ~ 

or equivalently, since e n R\ — > a, 

lim Rl F S JZ £ " hits R n A c ) = o«(0) = u (Ka) (0). 



n—too 

?2 



This completes the proof of ([35]), in the case a > 0. □ 
Proof of Lemma By replacing A with bA, we may and shall assume in this proof that 6 = 1. We 
thus have b n — > 1 as n — > oo. We first prove that 

liminf P r _i 1A (X £ " hits b n A c ) > P S JY^ hits A c ). (38) 

By Lemma[3j the events {y( K ) hits A c } and {y( K ) hits (^4) c } coincide a.s. We can then find a countable 
collection {ipi)i>\ of continuous functions with compact support contained in (A) c , such that 

{y<"> hits (A) c } = { sup (sup (Yt K \<pi)) > o), P So a.s.. 
i>\ v t>o 7 J 

Hence, if (tj)j>i is a sequence dense in [0, oo), we have 

P So (Y^ hits (A) c ) = lim tP,5o( sup f sup <Yf\^}) > 0). (39) 

However, Theorem [2] implies that, for every N > 1, 

liminf P r V ( sup ( sup {Xg, &)) > o) > P So ( sup ( sup (# , ^)) > o) . (40) 

Recall that b n — >■ 1, and note that the support of each function ^ is at a strictly positive distance 
from the set A. As a consequence, for every fixed N, the support of (p^ will be contained in b n (A) c for 
every i = 1, . . . , JV, as soon as n is large enough. Hence, for all large enough n, 

{ sup ( sup (X £n ,<pi)) > o) C hits b n A c }. 

l<i<N ^ l<j<N 3 ' > 
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Using this inclusion and then (|40p and (|39p . we immediately obtain (I38J) . 
We next turn to the more difficult upper bound 

limsup F,-i ]s (X £n hits b n A c ) < P 5o {Y (k) hits A c ). (41) 

We fix 5 > small enough so that the closed ball of radius 45 centered at is contained in A. As in 
Lemma |2j we let As be the connected component of the open set 

{x <E A: dist(x,^ c ) > 5} 

that contains 0. We denote the exit measure from As for the rescaled branching Brownian motion X Sn 
by £g. In other words, the measure ££ is equal to e n times the sum of the Dirac point masses at all 
points of OAs which are first exit points from As for one of the historical paths associated with X £n 
(these historical paths are defined in Subsection 12,31 for the branching Brownian motion Z 6n , and this 
definition is extended to X £ " by an obvious scaling transformation). 

Let $ be a continuous function on M, d such that 0<<I ) <l, < l ) = 0on ^3,5 and $ = 1 on AZg- Then, 
for every 77 > and p > 0, 

¥ le^]S ( X£n hitS b nA C ) =P [e -i ]5o (X e " hits b n A C , (£2, 1) < T/) 

+ P [£ -i ]5o (V" hits b n A c , (£%,!) > r,, J°°{X?,$)d8 < 

+ F [e - 1]s (x^ hits b n A c , (££,1) > r,, J™(X?,*)da > p\ (42) 

Let a n (rj), (3 n (j],p) and J n (r],p) be the three terms appearing in the right-hand side of (l42j) in this 
order. 

We first bound a n (rj). Provided n is sufficiently large, b n A c is contained in Ag, 2 and thus 

a n ( V ) < ^]S ( X£ " hits A h> > 1) < V)- 
Note that the times at which the historical paths of X £n exit As form a stopping line in the sense of 
[Ch9lj. We can thus apply the strong Markov property at a stopping line (Proposition 2.1 in [Ch91j. 
or more precisely a slight extension of the results in |Ch91| since our spatial motion is not standard 
Brownian motion, but Brownian motion killed inside T ro ) to obtain that 

¥ le^]s ( x£ " A c s/2 , (£^,1) < r,) =E [£ - 1]fio [l {(£ „ 1)<7?} P^(X £ « hits A c s/2 ) 
where £ V 1 is the point measure defined by £f(B) = e~ x £^ (e]l 2 B) for every Borel subset B of Wl d (££ 

— 1/2 

is in fact the exit measure of the unsealed process Z £n from e n As). From this equality, we see that 
a n (77) is bounded above by the probability for a branching Brownian motion (without killing) starting 
initially with less than r]e~ l particles, that one of the historical paths reaches a distance greater than 
5/(2yfe^) from its starting point. The estimate ([8j) now gives 

a n (v) < Vn 1 x n (^ (0) hits B(0, 5/(2^)) c ) < C'((d, v) ^. (43) 

Then, we have 



Recall that $ = 1 on A 2 g and in particular $ = 1 on B(x,5) for every x S dAs- We use the strong 
Markov property at the same stopping line as in the previous argument, together with a simple coupling 
argument, to write that 



Pn(T),p) < F [vE - 1]So ^j^(X^,l^ m )ds < , 
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where X £n is defined in terms of a branching Brownian motion Z £n in the same way as X £n was 
defined from Z £n . This branching Brownian motion Z £n has the same offspring distribution as Z £n , 
but particles are now killed at rate e n homogeneously over M. d . Furthermore, Z £n also starts from Mo 
under the probability measure Pfc<5 - 

By Proposition [2j the law of (X^ n )t>o under Pr^-i]^ converges as n — > oo to the law of Y^ 1 ' under 

P^Sq (in the notation of Subsection 12, 4p , Noting that, for every fixed s > 0, Yg a.s. does not charge 
the boundary of the ball B(0, 5), it follows that 

lirn^supP^- 1]5o Q"(^",l S(0i5) )d S < < P^y^iYW^^ds < p^j =: &ofa.P)- 

The continuity of sample paths of Y^ 1 ) ensures that /3oo(f? ; p) — ► as p — > 0, for every fixed r\ > 0. 
For the term 'yni'HiP)-, we simply use the bound 

7n(r?, p) < P [e -i ]5o ( j^(X £ s -,^)ds > p^j . 
This bound and the weak convergence of Theorem [2] imply that 

hm sup 7 „ < P So ( r(Y} K \<S>)ds >p)< P 5o (Y ( ^ hits A c 3d ), 

n— >oo ^ JO 

since <3? = on ^3,5. (To justify the first inequality in the last display, we also use the fact that the 
extinction times of X £n under P, -1,^ are stochastically bounded, which follows from a standard result 
in the case without killing.) 

To complete the argument, fix -d > 0. By Lemma [2] (ii), we can choose 5 > sufficiently small so 
that 

P So (Y^ hits A C 3S ) < P 5o (Y^ hits A c ) + |. 
From 043 p . we can then choose 77 > sufficiently small so that for all large n, 

Finally we choose p > such that /3oo(r/, p) < ^. From (|42p and the previous estimates, we obtain 
limsup P fe -i la (X £n hits b n A c ) < Ps {Y iK) hits ,4 C ) + 0, 

and since $ was arbitrary this completes the proof of (141 [) and Lemma [3 □ 
We still have to discuss the case a = in 035 p . So, let us consider two sequences {e n )n>i and 
(R n ) n >i such that e n i?^ — > 0. Let oq > and e' n = e n V (aoR~ 2 ). Since e n < e' n , we have 



liminf i£ F 5 JZ £n hits i? n A c ) > lim inf R 2 n ¥'$,,(Z £n hits i? n A c ) = U( Kan) (0), 

ri— >oo n— >oo 



by the case a > 0. By Lemma [2] (i), it( Kao )(0) can be made arbitrarily close to it(o)(0) when 00 is small, 
and so 

liminf R 2 n W 8 ,AZ £n hits R n A c ) > u (0) (0). 



To obtain the corresponding upper bound, a similar coupling argument shows that it suffices to consider 
the case when e n = for every n, that is when there is no killing inside the obstacles. Hence, consider 
the branching Brownian motion Z° = (the notation is even more legitimate since Z m,a does not 

depend on w). For every p > 0, define a rescaled version of Z° by setting 
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By Proposition [2j the law of (X t P )t>o under Pr^-ii^ converges to the law of under P§ as p tends 
to 0. Set p n = Rn 2 , and note that 

{Z° hits R n A c } = {X ipn) hits A c }. (44) 

A simplified version of the arguments of the proof of Lemma [7] shows that 

limsupPr -i ]s (X ipn) hits A c ) < P So {Y (0) hits A c ) = 1 - exp(-u (0) (0)). 

Arguing as in the first part of the proof of the theorem and using (j44|) yields 

limsupi?^P 5o (Z £n hits R n A c ) < u (0) (0), 

n— >oo 

which completes the proof of Theorem [1] □ 
Proof of Corollary [H For every r > 0, 

V neSo (V^R™' £ <r)=(l- ¥ So (Z^ £ hits B(0,e- 1 /2 r ) c ))" B . 
However, Theorem [1] shows that, P(cfaj) a.s., 

2 

n e F SQ (Z m > £ hits B(0, e- l ' 2 r) c ) = ^ x (?-P So (Z m > £ hits S(0, e _1/2 r) c )J 

converges to ^ n° Kj , 2 )(0) as e — > 0. The desired result follows. □ 
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